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1.  GENERAL  INTRODUCTION 

t 

x — The  interatomic  forces  in  a  crystal  are  usually  strongly  depen¬ 
dent  on  the  Interatomic  spacing  and  consequently  the  lattice  potential 
energy  for  a  crystal  can  generally  be  written  as  a  power  series  in  the 
displacements  of  the  atoms  from  their  equilibrium  positions.  If  such  an 
expansion  is  terminated  at  the  quadratic  terms,  then  this  constitutes  the 
so-called  harmonic  approximation.  In  the  harmonic  approximation  the  lattice 
vibrations  are  true  normal  modes  such  that  if  energy  is  channeled  uniquely 
into  any  one  lattice  mode,  then  it  will  remain  undissipated  in  that  mode. 

The  harmonic  approximation  would  theretore  predict,  for  instance,  that  the 
scattering  cross  section  for  neutrons  by  a  Bravais  lattice  would  consist  of 
a  set  of  6  function  peaks  or,  analogously,  that  the  spectral  profiles  of 
lattice  vibrations  would  be  a  set  of  undamped  temperature  independent  reso¬ 
nances.  These  and  other  such  predictions  are,  of  course,  in  marked  discord 
with  experimental  results.  It  is  therefore  clear  that  the  anharmonic  terms 
In  a  lattice  potential  must  be  considered  if  a  full  understanding  of  many 
of  the  physical  properties  of  solids  is  to  be  achieved. 

The  inclusion  of  the  anharmonic  terms  m  a  tQt  ice  potential  has 
several  consequences  as  far  as  phonons  are  concerned.  Firstly,  the  phenom¬ 
enon  of  thermal  expansion  is  now  allowed  and  this  leads  to  so-called  thermal 
strain  shifts  of  the  phonon  energies  away  from  their  harmonic  values. 
Secondly,  interactions  between  the  normal  modes  can  now  occur  and  this  opens 
up  channels  for  the  decay  of  phonons  which,  in  turn,  causes  additional  so- 
called  anharmonic  self-energy  shifts  of  phonon  energies  away  from  their 

(1  2) 

harmonic  values  and  the  appearance  of  finite  lifetimes  for  the  phonons.  ’ 
Thirdly,  it  opens  up  the  possibility  of  interactions  between  certain  phonons 
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(3  4) 

and  fluctuations  in  the  phonon  density  distribution,  *  one  of  the  ramif¬ 
ications  of  which  (under  favorable  conditions)  can  be  the  appearance  of  a 
relaxing  self-energy  component  for  a  mode,  and  it  is  this  latter  which  is 
predicted  to  lead  to  low-frequency  structure  in  the  response  function  for 
the  mode,  which  is  now  generally  referred  to  as  a  central  component. 

Before  proceeding  further,  it  is  useful  to  stress  certain  dis¬ 
tinguishing  features  and  certain  common  characteristics  of  these  components. 
The  thermal  strain  component  for  a  mode  is  a  frequency  independent  quantity, 
wheress  the  anharmonic  self-energy  and  relaxing  self-energy  components  cte 
frequency  dependent.  The  two  latter,  however,  are  differentiated  in  their 

frequency  dependence  because  the  main  structure  in  the  anharmonic  self-energy 
12 

occurs  around  10  Hz  whereas  the  central  component  arising  from  the  relaxing 

6  9 

self-energy  occurs  at  much  lower  frequencies,  probably  in  the  10  -10  Hz 
range.  The  ~entr3l  feature  which  is  common  tc  all  these  components  is  that 
they  each  can  contribute  to  the  temperature  dependence  of  a  mode,  somet'mes 
fn  a  dominant  and  Important  way.  In  weakly  anhaimonic  systems,  like  the 
alkali  halides  for  Instance,  the  temperature  dependence  of  all  modes  is 
generally  determined  by  the  competition  or  enhancement  between  the  thermal 
strain  and  the  anharmonic  self -energies ;  for  such  materials  the  magnitude 
of  both  components  for  a  mode  is  generally  small  compared  to  the  harmonic 
energy  of  the  mode,  so  that  the  resulting  temperature  dependence  is  also 
small.  The  effect  of  these  two  contributions,  however,  can  be  separated 
because  the  probability  for  phonon  interaction  is  determined  by  the  thermal 
population  factors  with  the  result  that  a  phonon  anharmonic  self-energy  will 
contribute  a  temperature  dependence  even  under  isochoric  conditions,  which 
is  in  contrast  to  that  arising  from  the  thermal  strain  component.  For 
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certain  modes  in  more  strongly  anharmonic  systems,  such  as  the  soft  mode  in 
displacive  ferroe lec tries ,  this  marginal  imbalance  between  the  magnitude  of 
the  thermal  strain  and  anharmonic  self -energy  is  broken,  with  the  latter 
becoming  so  large  as  to  strongly  dominate  both  the  thermal  strain  and  the 
mode  harmonic  energy,  which  is  imaginary,  so  that  such  systems  are  stabilized 
in  their  paraelectric  phase  by  these  huge  anharmonic  self-energy  contribu¬ 
tions  to  the  soft  mode.  In  certain  other  strongly  anharmonic  systems,  like 
piezoelectric  hydrogen-bonded  ferroelectrics ,  the  relaxing  self-energy  can 
also  wontribuce  importantly  to  the  temperature  dependence  of  the  soft-mode 
in  that  it  can  cause  the  mode  to  condense  out  below  the  characteristic 
clamped -Curie  temperature.  In  addition,  the  relaxing  self-energy  gives  rise 
to  a  temperature  dependence  of  the  central  component  which  is  most  pronounced 
as  the  clanq>ed-Curie  temperature  is  approached. 

These  three  components,  the  thermal  strain,  the  complex  anharmonic 
self-energy  and  the  relaxing  self-energy  provide  the  general  focus  of  the 
scope  of  this  final  report.  This  focus  is  inexorably  linked  with  the 

ultimate  goal  of  achieving  a  satisfactory  understanding  and  characterization, 
via  theoretical  descriptions,  of  anharmonic  interactions,  both  weak  and 
strong,  and  the  role  that  these  pla/  in  systems  like  ferroelectric  materials. 
Ferroelectric  materials  are  of  great  interest,  of  course,  because  of  the 
enormous  diversity  and  inter-relation  of  their  physical  properties  which  may 
be  harnessed  to  provide  solutions  to  technological  problems  such  as  memory 
function  devices  and  Che  like. 

Section  II  of  this  report  is  concerned  with  the  determination  of  the 
thermal  strain  component  and,  more  importantly,  of  the  frequency  dependence 
of  the  anharmonic  self-energy  and  damping  function  of  the  q  «  0  transverse 
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optic  mode  in  simple  ionic  systems.  Section  III  deals  with  an  experimental 
determination  of  the  self-energy  contributions  associated  with  the  soft  mode 
in  displacive  ferroelectric  materials.  Section  IV  describes  optical  and 
dielectric  investigations  of  the  relaxing  self-energy  associated  with  the 
soft  mode  in  piezoelectric  hydrogen -bonded  ferroelectrics . 
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II.  WEAKLY  AN HARMONIC  SYSTEMS 
A.  Introduction 

In  describing  the  physical  properties  of  a  real  anharmonic  crystal 
it  is  usual  to  resort  to  the  use  of  perturbation  theory.  While  ordinary 
perturbation  theory  has  been  used  to  describe  anharmonic  interactions,^  ^ 
more  recent  discussions  have  appealed  to  the  techniques  of  quantum  field 
theory  in  whion  the  system  is  described  in  terms  of  the  Green's  functions 
or  propagators  for  the  system.  The  use  of  temperature-dependent  time- 
ordered  Green’s  function*  to  describe  anharmonic  interactions  between  normal 
modes  of  vibration  has  been  described  by  Maradudin  and  Fein^^  and  Cowiey.^^ 
These  authors  have  shown  that  the  dielectric  and  scattering  properties  (which 
are  of  interest  here)  of  an  anharmonic  crystal  are  dependent  on  the  Fourier 
transforms  of  certain  time  correlation  functions,  the  simplest  of  which  is 
the  one-phonun  Green's  function  defined  by 


G(q,jj',t)  -  <TA(q, j,t)A*(q, j ',0)> 


vD 


where  the  phonon  operator  A(q,j)  is  defined  in  terras  of  the  sum  of  a  creation 
and  destruction  operator,  T  is  the  Dyson  time-ordering  parameter,  and  the 
triangular  parentheses  represent  thermal  averaging.  These  Green's  functions 
are  periodic  in  the  complex  time  direction  and  can  be  expanded  in  a  Fourier 
series  in  that  direction.  The  coefficients  of  this  series  are 


-ifift 


G(q,jj',n)  -  ^  J  G(q,jj',t)  eiC*  dt  , 


(2) 


where  £  -  1/kT,  k  being  Boltzmann's  constant,  and  Q  *=  2nzi/Sft,  x  being  an 
Integer.  The  physical  properties  of  the  crystal  can  be  obtained  from  these 
coefficients,  analytically  continued  over  the  whole  of  the  complex  SI  plane . 
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The  coefficients  are  obtained  from  the  use  of  diagrammatic  perturbation 
theory  and  the  Dyson  equation  for  the  Green's  functions  obtained  from  these 
diagrams  is  the  matrix  equation 

£  <{[u>h<q,j)]2  "  n2}  6 . , «  +  2<Uh<q,j)D(q,jj'.n)> 

i' 

6..*  2ujh(q , j) 

XGCq.jY.fl)  -  ~ -  ,  (3) 

where  DCq.jj'.O)  is  the  anharmonic  self-energy  of  the  phonons.  In  the  har¬ 
monic  approximation  the  «elf-energy  D(q,jj,,n)  is  zero  and  the  left-hand 
side  of  equation  (j)  then  plays  the  role  of  the  dynamical  matrix  within  the 
harmonic  approximation.  In  the  anharmonic  crystal  the  self-energy  is  not 
zero  and  the  matrix  on  the  left-hand  side  of  equation  (3)  couples  phonons 
from  the  dispersion  branches  j,  j',  and  j'  with  the  same  wave  vector.  This 
coupling  will  occur  whenever  the  modes  transform  according  to  the  same  irre¬ 
ducible  representation  of  the  space  group  of  the  crystals.  The  matrix 
has  Hermitian  and  anti-Hermitian  parts,  the  Hermitian  parts 
giving  rise  to  a  shift  in  the  normal  mode  frequencies,  and  the  anti-Hermitian 
parts  giving  rise  to  a  now  finite  lifetime  for  the  phonon  state.  Formally 
the  anharmonic  contribution  to  the  Hermitian  part  of  the  self-energies  of 
the  normal  mooes  can  be  included  by  renormalizing  the  frequencies  and  eigen¬ 
vectors  of  the  normal  modes.  If  the  off-diagonal  Hermitian  terms  in  the 
matrix  equation  are  neglected,  then  the  Green's  function  for  the  anharmonic 

crystal  becomes  similar  to  that  of  the  harmonic  crystal  if  [uj  (q,j)]  is 

2 

replaced  by  tu>(q,j)]  such  that  1 

C<“(q»j)]2  "  [u>h(q.j)]2  +  2u>h (q,j)D(q, jj.fi)  ,  (4) 
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where 

D(q,Jj,0)  -  Mq.Jj.n)  -  ir(q,jj,n)  •  (5) 

Detailed  expressions  for  the  self-energy  D(q,jj',Q)  have  been 
evaluated  by  a  number  of  authors  / 1,2 ,8,9)  To  second  order  in  both  cubic 

and  quartic  anharmonic! ty,  the  real  part  of  the  total  anharmonic  self-energy 
nay  be  written 

Mq.jj'.O)  -  AE(q,jj')  +  AA(q,jj\C0  ,  (*> 

where  AE(q,jj')  Is  the  thermal-strain  contribution  given  by 

AE(q,Jj')  «  |  *  Vqj;  -qj')  u^  ,  (7) 

and  &T(q,jj',fl)  is  the  anharmonic  self-energy  contribution  given  by 

iA(q,Jj'.fi)  -  AA<q,Jj'>  +  AA(q,jj',n)  ,  (6) 

where  A(q.jj')  is  a  frequency- independent  contribution  and  A  (q»jj  >^)  13 
a  frequency-dependeut  contribution  to  A  (qfjj  given  by 

AA(q,jj')  **  A^(q.Jj')  +  A^(q.jj')  * 

with 

A<4'(q,jj')  *  T  S  V(qj;  qj';  qljli  ”qlJl)<2rll  ^  l)  ’  (1°' 

qi  h 

A(8)(q,jj')  -  ^  ^  V(qj;  qj';  -q^S  -q^) 

*  ,Sl^l  ^2  -3^3 

+  n2  +  1  nx  -  nA 

XV(q.J.;  -<hV  *hV  "q3j3\  h  h  h 

~  u>i  +  u>2 

X  (2nj  +  1)  ,  (11) 


8- 


and 


-  A<6)(q,jj',n)  +  A<8)(q,Jj',n)  , 


(12) 


with 


d(6)(q.jj',n> 


-t!  2 ,  f  lv<S^  Si h>  -aiVl2  R(n) 

*  Sl^l  J2 


(13) 


A(8)(q,jj',n)  »  -Q  r  E  Z  |vfqj;  q^;  q2J2;  q3j3)|2S(n),  (14) 
*  Si-5!  S2J2  S3-3 


where 


R(0) 


+  n2  +  1  n1  +  ^  +  1  2(11!  “  n2^ 


U  +  0), h  + 


h  _  h  h  _  h  . 

j  +  tt>2  n  -  ^  -  o»2  n  -  oji  +  u2 


£>  • 


(15) 


SW>  -  [[(nl  +  1)(n2  +  1)<n3  +  l)  ■  nln2n3](n  ““7~h  "  h 

W  +  cu,  +  ux  +  ux 


1  2  3 


h  .  ..  h  _  h)+  3rr‘l(n2  +  1)(n3  + 


nu  li 

’  *1  "  ®2  "  ®3 


-  (nj  +  h  h  h 


«/  "Vv  i  *  .  •»  »  n  a  * 

Ti  -  U/j  +  u)2  +  U)^  0  + 


h  h 

iBj  -  -  aij 


x>]- 


(16) 


The  imaginary  part  of  the  self-energy  is  given  by 


r(q,)j',n)  «  r(6)(q,ij',n)  +  rw(q,jj',n)  , 


(8), 


(17) 


where 


SA(q,jj',n) 


!  p  r(q,jj',n) 
*  J..  (n-n;)p 


dfl' 


(18) 


In  these  equations  the  n's  are  phonon  population  numbers,  the  V  coefficients 
are  the  Fourier-transformed  anharmonie  force  constants,  and  the  superscripts 
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on  the  A  and  T  label  that  part  of  the  quantity  coming  from  the  contributions 
to  that  order  in  parameter  when  the  Hamiltonian  is  of  the  form 

H  *=  Hq  +  7]3H3  +  •  (19) 

In  this  section  of  the  report,  we  are  concerned  <?ith  the  theoreti¬ 
cal  and  experimental  determination  of  the  quantities  defined  by  equations 
(7),  (8)  and  (17)  determined  at  q  0  and  with  j  •=  j'  =  t,  chat  is  for  zone- 
center  transverse  optic  phonons  in  simple  cubic  ionic  systems  like  the 
alkali  halides. 


B.  Numerical  Evaluations  of  Self-Energy  Components 

As  can  be  seen  from  the  previous  section,  the  self-ene'gy  com¬ 
ponents  of  interest  here  involve  complex  lattice  dynamical  summations  which 
require  significant  computational  times  even  on  today’s  high-speed  computers. 
As  a  result,  reports  of  such  calculation*  havp  been  somewhat  limited  and  to 
date  have  been  confined  to  a  few  simple  ionic  materials. 

Such  numerical  estimates  of  the  different  self-energy  components 
can  be  made  providing  the  harmonic  frequencies  U)  (^,  j)  are  known  throughout 
the  Brillouin  zone  and  provided  the  J  coefficients  can  be  evaluated.  For 
simple  cubic  materials  like  ‘he  alkali  halides  reasonable  lattice-dynamical 
calculations  of  the  dispersion  curves  throughout  the  Brillouin  zone  can  be 
made  thus  generating  the  eigenvectors  and  eigenvalues  necessary  for  the 
evaluation  of  the  lattice  sums  contained  in  equations  (7)-(18).  Alterna¬ 
tively,  these  can  be  obtained  by  fitting  measured  iow-temperaiure  dispersion 
curves  to  particular  lattice-dynamical  models.  In  order  to  evaluate  the  V 
coefficients,  certain  assumptions  must  be  made  concerning  the  nature  of  the 
forces  that  interact  between  the  ions  in  the  lattice.  The  procedure  gener- 


t 
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ally  chosen  is  to  assume  a  short-range  central-force  repulsive  potential 
and  to  ignore  any  anharmonicity  in  the  Coulomb  contributions  to  the  lattice 
potential  because  of  the  expected  dominance  of  the  former  contribution  (the 
anharmonicity  of  the  Coulomb  interactions  can  be  included  if  necessary, 
though) . 

The  specifics  of  the  various  lattice  dynamical  calculations  that 
are  used  to  generate  the  harmonic  frequencies,  j),  and  of  the  procedures 

to  evaluate  the  V  coefficients  are  lengthy  and,  since  they  have  been  dis¬ 
cussed  in  the  literature , ^ ^  we  will  not  detail  them  here. 

E  A 

This  group  has  made  numerical  calculations  of  A  (0,t),  A  (0,t,U) 
and  r(0,t,n)  and  their  temperature  dependence  for  all  of  the  alkali  halides. 
In  making  our  calculations  of  the  principal  parts  and  6  functions  occurring 
in  equations  (7)-(18),  we  used  the  representation  method  suggested  by 
Maradudin  and  Fein,^  which  involves  the  approximation 
'  C(x)  3'1  '  -  (*  +  ie)'1  ,  (20) 

where  e  has  a  small  but  finite  value.  In  our  calculations^^'*^  we  used  a 
mesh  of  8000  points  in  the  Brillouin  zone  for  r^(0,t,Q)  and  A^(0,t,Q), 

✓  g\ 

but  only  a  tnesh  of  1000  points  for  the  lengthier  calculations  of  VK  0,tff5) 

/  o\ 

and  Av  (0,t,Q).  In  general  we  used  a  value  of  e  -  0.01  u  ,  vhere  tu  is  the 

L  L 

value  of  the  highest  frequency  in  the  Brillouin  zone. 

C.  Experimental  Determinations  of  Self-Energy  Components 
(1)  Introduction 

Under  normal  conditions,  experimental  measurements  do  not  lead 
straightforwardly  to  a  separation  of  the  thermal  strain  and  anharmonic  self¬ 
energy  components.  However  ,  as  Lowndes  and  Martin^7^  first  pointed  out, 


% 
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they  can  be  separated  out  by  a  suitable  combination  of  the  constraint  of 

high  pressure  and  variable  temperature  because  the  anharmonic  self-energy 

components  will  contribute  a  temperature  dependence  to  a  phonon  energy  that 

occurs  even  when  a  crystal  is  maintained  under  isochoric  conditions.  Lowndes 

and  Martin^^  used  this  approach  via  high-pressure/variable  temperature  low- 

3 

frequency  10  Hz)  dielectric  constant  measurements  to  make  first  estimates 

E 

of  the  thermal  strain,  A  (0,t),  and  low-frequency  anharmonic  self-energies, 

A 

A  (0,t,0),  of  the  ^  «  0  transverse  optic  phonons  in  simple  ionic  solids.  In 
the  course  of  the  work  covered  by  this  final  report,  this  group  has^^’^^ 
developed  high-pressure  far  infrared  techniques  in  order  that  a  similar 
approach  via  high-pressure/variable  temperature  far  infrared  measurements 

could  be  made  to  make  separate  estimates  of  the  thermal  strain  component, 

E  12 

A  (0,t),  and  first  estimates  of  the  far  infrared  (~  10  Hz)  anharmonic 

A 

self-energy  components,  A  (0,t,«a  ),  of  the  q  as  0  transverse  optic  phonons  in 
the  same  simple  ionic  compounds. 

But,  as  equations  (8)  and  (17)  show,  the  anharmonic  self-energy 
components  are  frequency  dependent  and  their  full  frequency  dependence 
really  needs  to  be  determined  in  order  to  more  fully  test  the  current  anhar¬ 
monic  LneorJes.  To  achieve  this,  one  needs  in  principle  to  be  able  to 
measure  the  frequency  dependence  of  both  the  real  and  imaginary  parts  of  the 
mode  response  function  via  some  convenient  physical  property.  Unfortunately, 
most  measurements  do  not  lead  directly  to  the  form  of  both  the  real  and 
imaginary  part  of  the  system  response  (conventional  power  far  infrared  spec¬ 
troscopy,  for  example,  leads  directly  to  only  the  reflection  amplitude, 
whereas  the  phase  angle  needs  also  to  be  determined  in  order  to  evaluate  the 


real  and  imaginary  pares  of  Che  dielectric  response;  similarly,  inelastic 
light-scattering  techniques  lead  directly  to  only  the  imaginary  part  of  the 
susceptibility).  Of  course,  in  principle,  such  measurements  can  be  used  to 
calculate  the  other  component  via  an  appeal  to  causality,  such  as  when  the 
Kramers -Kronig  relations  are  used,  but  such  procedures  have  only  limited 
accuracy.  One  measuring  technique  which  overcomes  these  problems  is  the 
comparatively  new  and  novel  technique  of  dispersive  Fourier  transform  spec¬ 
troscopy^^  (DFTS),  in  which  a  sample  is  placed  in  one  arm  of  a  Michelso 
interferometer  rather  than  outside  it  as  in  conventional  spectroscopy.  DFTS 
leads  directly  to  a  simultaneous  determination  of  the  system  reflection 
amplitude  and  phase  angle,  and  hence  of  the  real  and  imaginary  parts  of  the 

dielectric  response.  This  in  turn  can  then  be  used  to  make  determinations 

A 

of  the  frequency  dependence  of  A  (0,t,Q)  and  r(0,t,Q).  During  the  work 
covered  by  this  final  report,  this  research  group,  in  collaboration  with 
Dr.  T.  J.  Parker  of  the  University  of  London,  has  developed  and  used  this 
technique  to  lead  to  a  determination  of  the  full  frequency  dependence  of 
AA  (0,t,Cl)  and  r(0,t,Q)  for  a  number  of  simple  ionic  materials. 
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(2)  Determination  of  AE(0.t) 

The  current  theories  of  anharmonicity  lead  to  the  conclusion  that 
the  quasi-normal  frequency  of  a  phonon  measured  in  any  resonance  experiment 
is  given  by 

i«r  0^0,t)  53  +  2u)No,t>iAA‘(0,t,u;r)  +  AE(0,t)}  =  u>t  ,  (21) 


where  the  subscripts  T  end  0  denote  the  operative  temperature  and  pressure 

f\/)  £ 

of  the  system,  respectively.  It  is  readily  shown  that  &  ( 0,t )  is 
determined  by 


A (0,t)  « 


K.p(0>t)2  "  ^r,o(°»t>2] 


2cuoo^0,t^ 


(22) 


where  uJq  g(0,t)  ^'as  be,Jn  taken  as  a  reasonable  estimate  of  u)h(0,t),  and 

aij,  p(G,t)  is  the  lattice  vibration  frequency  measured  at  a  temperature  T 

and  at  a  pressure  P  which  is  such  as  to  reduce  the  volume  of  the  crystal 

to  that  which  it  has  at  0  K  and  zero  pressure.  The  right-hand  side  of 

equation  (22)  contains  all  experimentally  measurable  parameters  thus 

£ 

enabling  an  experimental  determination  of  A  (0,t). 

E 

A  separate  determination  of  &  ' 0 ,t )  can  be  made  as  follows.  At 
n  «  0,  r(0,t,0)  =  0,  and  the  lattice  contribution  to  the  static  dielectric 
response  can  be  written  as 


eL(0)  „ _ SW2 _  .  (2: 

tuh(0,t)2  +  2mh(0,t)UA(0,t,0)  +  &E(0,t)] 

2 

Providing  the  oscillator  strength  S(0)  is  reasonably  temperature  indepen¬ 
dent,  it  is  readily  shown^^  that  equation  (23)  leads  to 
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AE< 0,t)  *  - 


o>T  O(0,e2  AeL(0)]T 

U0,0(0  cJi07o) 


(24) 


where  the  subscripts  on  u  and  e  again  refer  to  temperature  and  pressure, 

respectively.  AeL(0)]T  is  a  measure  of  the  change  in  e^(0)  on  raising  the 

pressure  isothermally  at  a  temperature  T  so  as  to  restore  the  volume  of  the 

crystal  to  that  which  it  has  at  0  K  and  zero  pressure.  The  approximation 

sign  used  in  this  equation  stems  from  two  sources.  Firstly,  as  discussed 

above  we  have  used  values  of  the  lattice  vibration  frequency  determined  at 

close  to  T  -  0  K  and  at  zero  pressure,  w  (0,t),  to  approximate  uh(0,t). 

00 

2 

Secondly,  it  is  assumed  that  S(0)  is  reasonably  teropeiature  independent; 
justification  for  this  for  temperatures  up  to  500  K  nas  been  given  by  Lowndes 
and  Martin.  The  right  hand  side  of  equation  (24)  contains  all  experi¬ 

mentally  measurable  parameters  and  hence  can  lead  to  a  second  determination 


of  Ac‘(0,t) . 

*** 

£ 

In  our  work,  we  have  determined  A  (0,t)  from  both  equations 

(22)  and  (24)  and  have  generally  found  good  agreement  between  the  two 
approaches. 

E 

The  determination  of  A  (0,t)  by  these  two  separate  procedures  there¬ 
fore  requires  the  measurement  of  the  isobanc  (at  1  bar)  temperature  dependence 

of  w  (0,t)  and  the  corresponding  high-pres&ure  and  temperature  dependence  of 

X  ,u 

£ 

u>  (0,t)  to  determine  A  (0,t)  from  equation  (22),  and  the  additional  measure- 
T*P 

ment  of  the  isobaric  (1  bar)  temperature  dependence  of  the  dielectric  response, 

th  (0),  and  the  corresponding  appropriate  high  pressure  and  temperature  depen- 
x  »U 

L  E 

dence  of  e  (0)  to  determine  A  (0,t)  from  equation  (24).  The  lattice  vibration 


frequency  for  the  q  to  0  transverse  optic  mode  in  the  alkali  halides  is  most 
conveniently  studied  via  (conventional)  far  infrared  spectroscopic  trans¬ 
mission  measurements.  This  laboratory  has  developed  high-pressure  far 
infrared  techniques  sufficiently  so  as  to  facilitate  the  necessary  deter¬ 
minations  of  tt>TO(0,t)  and  uiT  p(0,t). 

(3)  Determination  of  AA  ( 0 ,  t  ,Ui  f)  and  r(0,t,UJt> 


At  T  =0K  and  at  zero  pressure ,  equation  21  becomes 


[<i>00(0,t)j2  =  [u)h(0, 

t)]2  +  2'Jih(0,t)&*(0,t,wc') 

(25) 

where  tt) t '  -  uQO(0,t) 

(26) 

while 

at  a  temperature  T  and  a  suitably  chosen  applied  pressure 

P  so  as  to 

maintain  the  crystal  volume  at  that 

found  for  the  crystal  at  OK 

and  zero 

pressure,  equation  21  becomes 

Cu)TjP(o,r)]2  =  Oh(0,t)]2 

+  >h(0,t)^(0,t,iut) 

(27) 

which 

leads  to 

Aj(0,t,«t)  -  d*(0,t,u>t') 

[tt)T  p(0,t)]2  -  L<»i00(0.t)]2 

(20) 

|! 

1 

6 

o 

o 

o 

rr 

’W 

1 

[*T.p(°,t)]?  -  [u)00(0,t)]2 

(29) 

s  -A'(T) 

HO) 

Equation  (30)  therefore  allows  an  experimental  determination  of 
the  change  in  the  anharmonic  self-energy  between  a  temperature  T  and  OK  via 

suitable  studies  of  the  temperature  and  pressure  dependence  oftti(0,t). 

A 

Aq( 0,t ,u>  ' )  is  not  necessarily  zero  because  anharmonic  interactions  can  per¬ 
sist  in  the  presence  of  zero  point  fluctuations.  However,  experimental  esti- 

A 

mates  of  Ag(0,t  can  be  made  as  follows.  For  temperatures  close  to  the  Debye 
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temperature  and  above,  it  can  be  shown  that  the  thermal  population  factors 

vary  linearly  with  T  and  also  that  the  leading  terms  in  both  the  cubic  and 

A 

quartic  anharmonic  contributions  to  A^CO.t.iu^)  are  linear  in  T.  If  such  a 
linearity  is  found  for  A'(T),  therefore,  extrapolations  of  these  linear 

A 

parts  back  to  T  *  0  will  yield  an  intercept  of  4Q(0,t,<Dt')  and  hence 

A 

A^.(0,t,U)t)  can  be  determined  from  equation  (29) 


(4)  The  Determination  of  AA(0.t,Q)  and  r(0,t,D) 

A  simple  cubic  material,  such  as  an  alkcli  halide,  has  a  single 

branch  of  transverse  optic  phonons  and  the  lattice  contribution  to  the 

dielectric  response  at  a  frequency  Q,  e^(fl) ,  is  proportional  to  the  Fourier 

transform  of  the  propagator  for  such  phonons  with  the  wave  vector  £  ss  0.  The 

Fourier  transform  of  the  propagator  leads  to  a  complex  lattice  contribution 

(1  2) 

to  the  dielectric  response  at  a  frequency  ft  of  the  form  ’ 


CL(Q)  =  _ U>h(0,t)2  S(fl)2 _ 

mh(Q,  r)2  +  2ujh(0,t)D(0,t,Q/  -  Q2 


(31) 


where  S(Q)  is  the  oscillator  strength  associated  with  thecas  0  transverse 

optic  mode  and  is  related  to  the  effective  dipole  moment  between  the  ions  in 

the  crystal.  The  mod?  anharmonic  self-energy  and  the  damping  can  be  expressed 

(21) 

in  terms  of  che  dielectric  response  of  the  system  as 


AA(o,t,n>  -  ^  {tuh(o,t)[s(n)V(n)  -  1]  +  - - )  *  aE(0,1)  >  (32) 

u>  (0,t) 


and 


<uh(Q.t)^"(Q)s(n)2 

2 


r(o,t,n) 


(33) 


where 


V<n) 


[e'(Q)  -  e(°»)] 


[e'(fi)  -  e(«)]2  +  e"(n)‘ 


(34) 


and 


TT'Cfl)  = 


[e'(n)  -  e(»)]  +  e"(n) 


2  ’ 


(35) 


where  e ' (Q)  and  e"(fl)  are  the  real  and  imaginary  part,  respectively,  of  the 

(total)  dielectric  response  of  the  system. 

Equations  32-35  reveal  that  experimental  determinations  of  the 
A 

frequency  dependence  of  A  (0,t,fl)  and  r(0,t,ft)  can  be  trade  if  valuer  of 
h  F 

(I)  (0,t),  S(Q),  T|’(ft)  and  T]"(fi),  and  A  (0,  t)  can  be  determined.  We  have  dis¬ 
cussed  the  determination  of  U)^(0,t)  and  A^(0,t)  in  section  C(2)  above. 

(22)  (21 

Bilz  et  al.  and  Cowley'  '  have  shown  that  the  oscillator 

strength  S (ii)  has  only  a  very  weak  frequency  dependence  compared  to  chat 

A  (is  23) 

arising  from  A  (0,t,fi)  and  r(0,t,fl).  Furthermore,  Lovmdes  and  Martin  ’ 

2 

have  sho/n  that  S(0)  has  only  a  very  weak  temperature  dependence.  At  T  =  0 

and  Q  =  0  the  dielectric  response  described  by  equation  (31)  reduces  to  a 

2  L  2 

good  approximation  to  S(0)  ca  e  (0)  and  so  values  of  S (0)  may  be  obtained 

to  a  good  approximation  from  values  of  e^(0)  which  have  previously  been 

(18) 

reported  by  this  author. 

Equations  (34)  and  (35)  reveal  that  the  frequency  dependence  of 

Tl'(n)  and  T)"(n)  are  determined  by  a  knowledge  of  a '  (Q)  and  e"(Q)  and 

(18  23) 

€(»)..  e(°°)  has  little  temperature  dependence.  *  We  have  used  disper¬ 

sive  Fourier  transform  reflection  spectroscopy  (DFTS)  to  determine  e '  (Q) 
and  €*'(0)  with  high  accuracy.  DFTS  will  allow  measurements  of  e ' (0)  and 
e"(Q)  to  be  made  simultaneously  to  an  accuracy  of  +  27.  throughout  the  far 
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infrared  region  (5-600  oq  V  It  is  to  be  noted  here  that  the  tempera¬ 
ture  measurements  of  e*(Q)  and  e*(fl)  need  to  be  determined  only  at  zero 

A 

pressure  in  order  to  determine  A  (0,t,Q)  and  r(0,t,Q). 

(5)  Experimental  Techniques 

A  variety  of  experimental  techniques  were  needed  to  measure 
the  various  parameters  required  in  the  determination  of  A  (0,t,n),  F(0>t>0) 
and  AE(0»t)  (see  equations  (32)-(35),  (22),  (24)).  aE(0,t)  was  deteiwired 

in  two  separate  ways:  through  (conventional  far  infrared  Fourier  transform 
spectroscopic  studies  of  the  phonon  frequency  and  through  studies  of  the 
low-frequency  dielectric  response,  both  made  at  high  pressures  and  variable 
temperatures.  The  determination  of  A  (0,t,fl)  and  r(0,t,ft)  required  the  use 
of  Dispersive  Fourier  Transform  Spectroscopy  to  measure  the  frequency 
dependence  of  the  complex  dielectric  response. 

(a) 

Cubic  ionic  solids  like  the  alkali  and  thallium  halides  have 
one  q  »  0  transverse  optic  mode  which  is  most  conveniently  studied  by  far 
infrared  spectroscopy.  The  frequencies  and  lifetimes  of  these  modes  are 
readily  determined  from  suitable  measurements  of  the  normal  incidence 
transmission  spectra  of  chin  films  of  the  material  under  investigation. 
Electromagnetic  analyses  of  thin-film  behavior  shows  that  the  minimum  ir.  the 
spectral  transmission  for  such  a  thin  film  occurs  at  u>(0,t)  s  providing 

I 

the  film  is  thin  compared  to  the  vacuum  wavelength  of  the  incident  radia- 
(24) 

Cion.  A  similar  analysis  reveals  that  2r(0,t,u>t)  is  determined  by  the 


(Conventional)  Fourier  Transform  Spectroscopy 
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spectral  damping  measured  at  a  transmission  T*,  where  T*  is  given  by 


1  +  T 


with  Tq  being  the  transmission  at  U)(0,t).  The  use  of  thin  film  transmission 
studies,  therefore,  allows  a  direct  determination  of  m(0,t)  and  2F(0,t,u)t) 
without  lengthy  analysis  of  the  data,  as  would  be  required  using  reflection 
data  from  bulk  samples. 

Values  of  uj^,  ^(O.t)  and  U)^  p(0,t),  which  were  required  to  determine 

g 

A  (0,t)  (see  equation  (22)),  were  measured  by  transmission  studies  on  thin 
films  using  conventional  far  infrared  Fourier  Transform  Spectroscopy. 

Values  of  («T  Q(0,t)  were  also  obtained  from  the  far  infrared  DFTS  measurements 
(to  be  discussed  beloj)  made  in  reflection  from  bulk  single  crystals  and 
hence  these  allowed  a  check  on  these  measurements .  It  is  to  be  noted  that, 
because  of  severe  technical  problems,  the  DFTS  technique  cannot  as  yet  be 
used  with  samples  under  pressure  and  hence  no  check  values  of  uiTp(0,t)  were 


obtained . 


The  measurements  of  uj(0,t)  were  obtained  by  means  of  Fourier 


Transform  Spectroscopy  using  Michelsoa  interferometers.  A  germanium  boivmeter 
(Infrared  Laboratories ,  Tucson)  operating  at  2  K,  was  used  as  a  detector  in 
the  frequency  range  3-300  cm  ;  a  Golay  cell  was  used  at  higher  frequencies. 


(b)  High  Pressure  Fourier  Transform  Spectroscopy 


The  requirements  of  hydrostatic  pressure  containment  and  far- 
infrared  spectroscopy  are  somewhat  conflicting.  The  first  requires  small 
apertures  for  maximum  strength  while  the  second  requires  large  apertures  to 
provide  maximum  energy  throughput  from  the  comparatively  weak  far- infrared 
sources.  These  requirements  are  further  complicated  because  of  the  scarcity 
of  readily  available  window  materials  which  satisfy  the  competing  demands 


-20- 


of  strength  and  transmissivity  in  the  far  infrared.  As  a  result  c  r  these 

difficulties,  much  of  the  very  limited  far-infrared  high-pressure  research 

has  been  confined  to  date  to  work  achieved  with  the  opposed  diamond  anvil 

(25) 

system.  When  used  properly  this  cell  is  a  powerful  tool  in  the  very- 

high-pressure  domain,  but,  in  the  important  0-10  Khar  range  involved  in  the 
present  work,  reliable  far-infrared  measurements  are  difficult  to  perform 
with  this  cell.  This  is  primarily  because  of  the  difficulty  in  accurately 
determining  the  system  operating  pressure  but  also  because  of  the  inherent 
difficulty  of  achieving  true  hydrostatic  pressure  conditions  because  of 
residual  pressure  gradients  across  the  anvil  faces.  The  consequences  cf 
these  problems  have  been  that  it  Is  difficult  to  precisely  determine  the 
pressure  dependence  of  mode  eigenfrequencies ,  especially  at  the  low  pressures 
necessary  for  the  evaluation  of  u^,  p(0,t),  and  it  is  virtually  impossible 
to  reliably  measure  the  (small)  pressure  dependence  of  the  associated  line- 
widths  using  the  diamond  anvil  system. 

Accordingly,  this  group  has  developed  a  high-pressure  far-infrared 
cell  which  goes  some  way  to  solving  these  problems,  1  and  this  is  illus¬ 
trated  in  figure  1.  The  cell  body,  window  mounts,  and  reraining  closures 
were  machined  from  4340  alloy  steel  and  then  heat  treated  to  a  hardness  of 
KC45  giving  a  yield  strength  of  14  Kbars.  The  critical  bores  and  surfaces 
of  the  cell  and  its  components  were  ground  and  honed  to  size  after  the 
hardening  process.  The  cell  body  measures  17  X  14  X  10  cm  and  fully  assem¬ 
bled  the  cell  has  an  effective  speed  of  f2.4. 

The  window  mounts  used  in  the  cell  give  a  supported  to  unsupported 
area  ratio  of  about  3  with  an  unsupported  port  diameter  of  5  mm.  The  window 
high-pressure  seal  is  achieved  by  lapping  and  polishing  the  mating  surfaces 
of  the  window  and  its  mount  to  be  flat  and  parallel  to  better  than  300A. 
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The  windows  are  prevented  from  sliding  off  their  mounts  by  brass  retaining 
rings  attached  to  the  window  mounts.  The  window  mount  pressure  seal  is 
achieved  using. a  stainless-steel,  brass,  and  Teflon  packing  assembly. 

In  choosing  both  the  cell  windows  and  the  pressure  transmission 
medium,  care  must  be  used  to  avoid  materials  which  have  characteristic 
electromagnetic  resonances  in  the  spectral  region  under  investigation.  The 
choice  of  materials  available  for  both  window  and  pressure  transmission 
media  for  far-infrared  service  is  extremely  restricted  nowever.  For  windows 
we  have  used  cj-in. -thick  fused  quartz  for  the  spectral  range  below  125  cm  * 
and  4- in* "thick  silicon  for  the  spectral  range  125-400  cm  For  a  pressure 
transmission  medium  we  have  used  either  helium  or  argon  gas  since  these 
exhibit  no  characteristic  far-infrared  resonances  as  do  most  of  the  organic 
oils  and  molecular  gases  conventionally  used  in  high-pressure  research. 

Figure  2  shows  the  experimental  arrangement  for  a  high-pressure 
far-infrared  transtr;i.ssion  experiment.  The  high-pressure  cell  is  located  in 
an  isolated  evacuated  thick-walled  module  B  designed  to  protect  the  inter¬ 
ferometric  and  detector  modules  A  and  C,  respectively.  The  cell  window  axis 
is  deliberately  chosen  to  be  off  line  with  both  the  entrance  and  exit  ports 
to  module  B  to  avoid  any  possible  damage  to  the  interferometer  or  detector 
due  to  a  high-pressure  window  failure. 

The  hydrostatic  pressure  for  these  experiments  was  generated 
from  a  two-stage  gas  compressor.  The  pressure  was  measured  using  a  suitably 
aged  and  calibrated  t«anganin  cell.  With  this  high-pressure  spectroscopic 
arrangement  we  were  able  to  make  spectral  studies  in  the  far  infrared  under 
truly  hydrostatic  pressure  conditions  with  a  pressure  accuracy  of  better 
than  +  17.  whilst  at  the  same  time  achieving  signal-to-noise  values  of  better 


than  100. 
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(c)  Low  Frequency  Dielectric  Constant  Measurement 

The  low-frequency  dielectric  response,  e(0),  was  determined  via 

three -terminal  capacitance  measurements  made  in  the  frequency  range  1-100 

(17  23  27) 

kHz  with  a  measuring  accuracy  of  ±0.17..  ’  *  The  capacitance  measure¬ 

ments  were  recorded  on  a  General  Radio  model  1615A  transformer-ratio-arm 
bridge  used  in  conjunction  with  a  tuned  amplifier  and  null  detector.  The 
measuring  electrodes  were  evaporated  onto  the  single  crystal  samples  using 
gold  as  the  electrode  material. 


(d)  Dispersive  Fourier  Transform  Spectroscopy 

A 

In  order  to  determine  the  frequency  dependence  of  A  (0,t,0)  and 
r(0,t,n)  from  equations  (32)  and  (33),  it  is  necessary  to  precisely  deter¬ 
mine  the  dispersion  of  e'(Q)  and  €"(0). 

The  conventional  method  for  determining  e,(0)  and  t" (0)  is  to 
measure  the  power  reflectivity,  R(fl),  at  near  normal  incidence  over  as  wide 

a  spectral  range  as  possible,  and  then  to  obtain  the  phase  spectrum,  cp(0), 

(28) 

by  calculation  from  the  Kramers -Kronig  relation' 


9(0) 


Xn  r(f)') 

-,2  2 
o  -  n 


dn' 


(37) 


In  this  expression  r(Q)exp[icp(fl)]  is  the  complex  amplitude  reflection 
coefficient  and  r(Q)  =  /r(D)  .  From  this  e'(fi)  and  eff(Q)  can  then  be  obtained 
from  the  Fresnel  relations.  The  main  problem  with  this  procedure  is  that 
the  complete  dispersion  of  R(0)  must  be  measured,  but  this  is  often  not 
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possible  especially  if  there  is  significant  dispersion  at  very  lcw- 

frequencies  (such  as  for  soft  modes)  or  at  high  frequencies,  which 

therefore  necessitates  other  measuring  techniques  to  be  utilized  to 

complete  the  measurement  of  R(fi).  A  second  problem  is  that  cp(fi)  is 

not  well  determined  in  those  frequency  ranges  where  r(Q)  *♦  0  because 

of  the  inaccuracies  in  measuring  the  magnitude  of  r  >'  f!) . 

Conventional  Fourier  transform  spectroscopy,  as  described 

in  a  previous  section,  yields  the  power  spectrum  from  a  sample 

placed  externally  to  the  interferometer .  However,  if  the  sample  is 

inserted  in  one  arm  of  the  Michelson  interferometer,  then  it  is 

possible  to  perform  amplitude-phase  or  dispersive  Fourier  transform 
(19  20) 

spectroscopy.  ’  The  asymmetric  sample -inter ferogram  from  such 

spectroscopy  can  be  transformed  into  the  spectral  functions  of  the 

phase  and  amplitude  of  the  reflectance,  from  which  the  real  and 

imaginary  parts  of  the  complex  dielectric  response  of  the  system  can 

be  measured  directly,  and  without  having  to  measure  the  infinite 

spectrum.  Furthermore,  since  DFTS  measures  r(Cl)»  rather  than 
2 

r  (0)  *  R(Q)  as  i.i  conventional  spectroscopy,  it  has  an  automatic 
improvement  in  the  accuracy  of  €7(n)  and  t" (fl)  which  is  especially 
important  as  r(Q)  -*  0. 

The  apparatus  used  by  this  group  is  illustrated  schematically 
(29-31) 

in  figures  3  and  4.  The  fixed  mirror  of  the  Michelson  interferom¬ 


eter  is  replaced  by  an  optically  flat  reflector  whose  surface  is  divided 
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up  into  an  outer  metallized  ring  C  and  two  inner  semicircular  rings,  one 
of  which,  B,  is  metallized  and  the  other  is  the  sample.  By  an  arrangement 
of  opaque  screens  each  of  these  parts  may  in  turn  be  used  as  the  reflector. 
In  practice,  the  outer  ring  is  first  used  to  align  the  whole  reflector  so 
that  the  incident  beam  is  perpendicular  to  the  surface.  Since  it  has  a 
large  radius  it  provides  a  readily  reproducible  alignment.  Then  a  specimen 
interferogram  is  recorded  with  part  A  exposed  and  a  reference  interferogrun 
is  recorded  with  part  B  exposed.  The  inter ferograms  are  Fourier  trans¬ 
formed  and  the  educed  spectra  are  ratioed  to  give  the  complex  reflectivity. 
To  reduce  systematic  errors  caused  by  asymmetry  between  the  two  sides  of 
the  instrument  this  output  spectrum  is  re-ratioed  against  a  ratioed  cali¬ 
bration  spectrum  obtained  in  a  separate  experiment  with  both  parts  A  and  B 
metallized. 

Although  the  specimen  is  on  a  collimated  beam,  the  method  is 

2 

suitable  for  measurements  on  fairly  small  samples  of  area  of  1  cm  or 
larger,  which  is  satisfactory  for  all  the  materials  studied  in  this  research 


project . 


The  procedure  of  metallizing  part  of  the  sample  surface  for  use 


as  h  phase  reference  surface  therefore  allows  the  sample  temperature  tc  be 

changed  and  overcomes  the  problems  experienced  by  earlier  workers  using  DFTS 

(19  20  ^2  34) 

which  confined  their  measurements  to  room  temperature  only.  * 

The  group  at  Northeastern  University,  in  collaboration  with  Dr.  T.  J.  Parker 
of  the  University  of  London,  has  used  this  technique  to  determine  e/(n) 
and  e"(n)  for  several  alkali  halides  (NaF,  NaCf,  K Cl,  KBr ,  KI,  RbCI,  RbBr 

A 

and  Rbl)  in  the  range  100-300  K,  and  hence  to  determine  A  (0,t,f2)  and 
(21,35,36) 

no,t,n). 
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D.  Results  and  Discussion 

A  E 

a)  Theoretical  Calculations  of  A  (O.t.Q).  r(O.t.fi)  and  A  (0,t) 

Figure  5  shows  the  frequency  dependence  calculated  for  the 

different  contributions  to  A  (0,t,n)  and  r(0,t,P.)  for  NaBr  and  Figure  6 

A 

shows  the  calculated  frequency  dependence  of  A  (0,t,P)  for  NaBr  at  dif¬ 
ferent  temperatures.  Such  results  are  generally  typical  for  all  of  the 
alkali  halides  crystallizing  in  the  NaCi  structure. 

£ 

Tables  1-4  summarize  our  calculated  results  for  A  (0,t), 

A  A 

A  (0,t,0),  A  (0,t,u)t)  and  r(0,t,u>t)  and  their  component  contributions 
for  the  lithium,  sodium,  potassium  and  rubidium  halides,  respectively. 

A 

The  data  for  A  (0,t,o>t)  and  r(0,t,ujt>  and  their  respective  components  are 
given  for  the  calculated  value  of  u^. 

An  analysis  of  the  data  in  Tables  1-4  reveals  a  number  of  over¬ 
all  qualitative  features  which  should  be  commented  on: 

A 

(i)  The  frequency-independent  contributi an  A  (0,t)  is  always 

dominated  by  the  positive  contributions  from  the  first-order  quart ic  term 

A^^(0,t);  A^^(0,t)  is  always  less  than  107.  of  A  ^(0,t)  even  above  the 

characteristic  Debye  temperature. 

~A 

(ii)  A  (0,t,0)  is  always  dominated  by  the  negative  second- 

order  cubic  term  A^(0,t,0).  Although  A^(0,t,0)  is  generally  less  than 
~A 

10%  of  A  (0,t,0)  at  low  temperatures,  it  has  a  much  stronger  temperature 

~A 

dependence  and  can  be  as  much  as  307.  of  A  (0,t,0)  at  the  characteristic 
Debye  temperature. 


S 
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(iii)  AA(0,t,0)  is  always  positive  at  low  temperatures  because  of 
the  dominance  of  the  first-order  quartic  term  A^)(0,t),  but  this  usually 
becomes  negative  at  higher  temperatures  because  of  the  stronger  temperature 
dependence  of  the  negative  contributions  of  A^(0,t,0),  a(®)(0,l)  and 
A<8)(0,t,0). 

(iv)  AA(0,t,wt)  is  always  dominated  by  the  negative  second-order 
cubic  term  A^8) (0, t,a>t) .  At  low  temperatures  A^®^ CO, t,wt)  is  always  less 
than  10%  of  AA(0,t,u)t)  but  because  of  a  strong  temperature  dependence  it 
can  be  as  much  as  40%  of  AA(0,t,wt)  at  the  characteristic  Debye  temperature. 

(v)  AA(0,t,tot)  can  be  either  positive  or  negative  at  low  tempera¬ 
tures  depending  on  the  balance  between  A^^(0,t)  and  A^8)  (0,  t.Uj.) .  At 
higher  temperatures  AA(0,t,wt)  is  generally  negative,  however,  because  of 
the  stronger  tempeiature  dependence  of  the  negative  contributions  from 
A^6)  (0,t,ut)  ,  A^8'  (0,t.»)t)  ,  and  A^®)(0,t). 

(vi)  r(0,t,wt)  is  totally  determined  at  low  temperatures  by  the 
second-order  cubic  terms,  r ^ (0, t ,uit) .  However,  i (0,t,wt)  has  a  much 
stronger  temperature  dependence  than  (O.t.u^)  and  both  contributions  to 
r(0,t,ut)  are  generally  very  close  in  magnitude  at  the  characteristic  Debye 
temperature. 

b)  Experimental  Results  for  AA(0,t,0),  AA(0, t ,u>t ) ,  rCO.t.w^)  and  AE(0,t) 

Using  the  appropriate  measurements  of  the  temperature  and  pressure 
dependence  of  w(0,t)  and  e(0),  we  have  determined  experimental  values  of 
AE(0,t),  AA(0,t,0)  and  AA(0,t,ut)  via  the  procedures  outlined  earlier,  and 


» 


these  values  are  summarized  in  Tables  5-10  for  the  lithi  n,  sodium,  potassium, 
rubidium,  cesium,  and  thallium  and  silver  halides,  respectively.  Table  11 
summarizes  the  temperature  dependence  of  2F(0,t,u)t)  for  all  these  materials. 

A  number  of  qualitative  trends  emerge  from  these  results:  (i)  the 
A 

dQCO.t.w^)  are  generally  quite  small  and  are  no  more  than  a  few  percent  of 

A 

u00(0,t).  (ii)  The  A^O.t.w^)  for  the  silver  and  alkali  halides  crystalliz¬ 
ing  in  the  NaCl  structure  are  generally  negative  in  sign  at  all  temperatures 

A 

and  increase  in  magnitude  with  increasing  temperature-  (iii)  The  A^,(0,t,u)t) 
for  the  thallium  and  alkali  halides  crystallizing  in  the  CsCl  structure  may 
be  negative  or  positive  at  lower  temperatures  but  are  always  positive  at 
higher  temperatures.  (iv)  The  A^(0,t,uj^}  have  a  smaller  temperature  depen¬ 
dence  than  AE(0,t)  for  all  the  materials  studied  except  the  thallium  halides. 

Tables  5-8  also  compare  the  theoretical  calculations  and  experimental 
results  for  the  anharmonic  self-energies  of  the  experimental  lithium,  sodium, 
potassium  and  rubidium  halides.  The  results  have  allowed  separate  decermina- 
tions  of  AE(0,t)  to  be  made  at  ft  =  0  and  ft  =  u>t.  Ecuation  (7)  shows  AE(0,t) 
to  be  a  frequency- independent  quantity  and L the  experimental  results  generally 
support  this  within  the  limits  of  the  experimental  accuracy.  The  exceptions 
to  this  are  the  results  for  the  cesium  and  thallium  halides  at  higher  tem¬ 
peratures,  but  these  differences  may  indicate  the  importance  of  the  volume 
dependence  of  the  dipole  moments  associated  with  the  transverse-optic  modes 
in  these  materials  which  have  been  neglected  in  the  determination  of  AE(C,t) 
at  ft  *  0.  It  is  for  this  reason  that  the  results  for  AE(0,t)  determined  at 
ft  ■  a»t  are  considered  more  reliable,  and  a  comparison  of  these  with  the  theo¬ 
retical  results  reveals  a  good  agreement  within  the  accuracy  limits  for  all 
the  materials. 

The  results  for  Aj(0,t,0)  are  qualitatively  very  good  in  that  both 
experiment  and  theory  show  a  reversal  in  sign  as  the  temperature  is  raised, 
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and  furthermore  show  this  reversal  to  take  place  at  about  the  same  tempera¬ 
ture.  Quantitatively,  the  low-temperature  results  are  in  good  agreement  for 
all  the  materials  and  this  reasonable  agreement  is  maintained  for  the  potas¬ 
sium  and  rubidium  halides  at  higher  temperatures.  For  the  lithium  and 
sodium  halides,  however,  serious  discrepancies  arise  between  the  two  sets 
of  data  at  higher  temperatures,  with  the  calculated  results  showing  a  much 
stronger  temperature  dependence. 

The  agreement  between  theory  and  experiment  is  less  good  for 

A 

A^(0,t,ut) .  Qualitatively,  the  experimental  results  suggest  that  the  second- 

A 

order  anharmonic  contributions  always  dominate  so  that  AT(0,t,u)^)  <  0, 
whereas  the  calculated  results  reveal  that  this  is  not  always  so  at  the 
lower  temperatures.  Quantitatively,  although  the  results  for  many  of  the 
compounds  agree  within  the  error  limitations,  the  experimental  results  gen¬ 
erally  tend  to  vary  somewhat  faster  with  temperature  than  the  calculated 

theoretical  predictions.  The  preceding  comments  have  been  made  for  the 

A 

calculated  values  of  A^,(0,t,iot)  determined  at  the  calculated  value  of  u)t. 

A 

In  order  to  test  the  sensitivity  of  the  calculated  values  of  AT(0,t,u)t)  to 

A 

any  frequency  dispersion,  we  have  also  listed  values  of  A^,(0,t,mt)  calculated 
at  the  experimentally  determined  values  of  in  Tables  5-8.  With  the  excep¬ 
tion  of  the  higher  temperature  data  for  LiBr  and  Lil  (which  in  rny  case  must 


be  treated  with  some  caution  for  the  reasons  given  earlier), the  two  sets 

A 

of  calculated  values  for  AT(0,t,wt)  are  quite  close  and  generally  differ  at 

A 

most  by  no  more  than  +  152.  However,  the  values  of  A^(0,t,wt)  calculated  at 
(u^)^  are  generally  in  no  better  agreement  with  the  experimental  values 
than  the  A^(0,t,wt)  calculated  at  ^ut)ca^c' 

The  results  for  r(0,t,wt)  show  the  theoretical  predictions  at  low 
temperatures  to  usually  be  somewhat  smaller  than  the  experimental  values 
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and,  as  the  temperature  is  increased,  to  generally  have  a  somewhat  faster 
temperature  dependence  than  the  experimental  results.  The  calculated 
theoretical  results  at  the  higher  temperatures  are  generally  in  reasonable 
agreement  with  the  experimental  results  for  the  lithium  and  sodium  halides 
but  are  significantly  bigger  than  the  experimental  results  for  the  potassium 
and  rubidium  halides. 

The  experimental  studies  of  the  pressure  dependence  of  rfO.t.uij.) 
Suggest  thac  its  volume  dependence  is  small.  Equations  (17)  and  (18)  re¬ 
veal  that  such  a  volume  dependence  can  arise  only  from  that  of  the  anharmonic 
coupling  coefficients  or  crom  that  of  the  c.h(q,j).  We  have  attempted  to 
calculate  approximately  what  Eqs.  (17)  and  (io)  would  predict  for  the  volume 
dependence  of  r(0,t,uO  for  a  number  of  alkali  halides  by  allowing  all  the 
w^(d,j)  to  have  the  same  pressure  dependence  as  that  measured  for  m(0,t) 
and  by  using  values  of  rQ  and  8  suitably  corrected  for  the  effects  of  pressure. 
Such  assumptions  lead  to  the  conclusion  that  r(Q,t,wt)  will  always  decrease 
with  increasing  pressure  by  (1-2)%  per  kbar.  Although  the  absolute  magni¬ 
tude  of  this  calculated  pressure  dependence  for  r(0,t,wt)  is  in  agreement 
with  that  found  experimentally,  these  results  do  not  explain  the  small  in¬ 
creases  with  pressure  found  for  r(0,L,«  )  for  some  materials. 

The  reason  for  the  less  good  agreement  between  the  calculated  and 
experimental  anharmonic  self-energies  at  ft  =  u>  is  not  clear.  The  discrep¬ 
ancy  between  calculation  and  experiment  for  r(0,t,o>t)  at  low  temperatures 
could  simply  arise  from  an  underestimate  of  r^(0,t,o>t)  or  (0,t,u>t)  or 
from  some  sample  artifact  in  the  low-temperature  experimental  measurements: 
this  latter  would  seem  most  unlikely,  however,  in  view  of  the  equivalent 
results  that  we  obtain  from  single-crystal  reflection  measurements  and  thin- 
filtn  transmission  measurements.  The  generally  good  agreement  between  cal- 


culation  and  experiment  for  the  higl rr  temperature  r(0,t,wt)  for  the  sodium, 
potassium,  and  rubidium  halides  suggests,  but  does  not  prove,  that  the  cal¬ 
culated  form  of  AA(0,t,ut)  may  be  correct;  if  this  were  true  then  this  in 
turn  would  suggest  that  the  general  disagreement  between  calculation  and 
experiment  for  Aj(0,t,ut)  might  arise  from  the  calculation  of  AA(0,t)  (see 
Eq.  (9)). 

c )  Experimental  Results  for  AA(0,t,Q)  an d  r ( 0 , t , H ) 

Figures  7-1h  show  the  experimental  determinations  of  the  frequency 
dependence  of  A^(0,t,fl)  and  r(0,t,fl)  for  NaCl  and  RbCl  at  105K  and  iOCK.  In 
all  we  have  measured  the  frequency  dependence  of  the  self-energies,  as  a 
function  of  temperature,  for  NaF,  NaCl,  KC1,  KBr,  KI,  RbCl,  RbBr  and  Rbl. 

Our  results  indicate  that  there  is  a  reasonable  qualitative  agree¬ 
ment,  within  experimental  error,  between  the  experimental  and  theoretical 
values  of  the  self-energies  for  all  materials.  Quantitatively,  the  results 
are  in  very  good  agreement  for  the  rubidium  and  potassium  halides,  but  there 
are  marked  discrepancies  between  the  two  secs  of  data  for  NaF  and  NaCl,  as 
the  figures  show.  It  should  be  noted  that  these  measurements  are  in  good 
agreement  with  our  earlier  measurements  of  the  self-energies  at  specific 
frequencies,  i.e.  AA(0,t,(P,  AA(0,t,uO  and  r(0.t,'nt). 

d)  Conclusion 

In  assessing  the  obvious  qualitative  and  quantitative  agreement  for 
the  calculated  and  experimental  determinations  of  these  anharmonic  self¬ 
energies,  it  is  useful  to  place  in  perspective  the  errors  associated  with 
each  kind  of  determination. 

Given  the  approximations  contained  within  the  anharmonic  theory,  the 
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errors  associated  with  the  calculated  values  of  the  anharmonic  self-energies 
stem  primarily  from  those  involved  with  the  ou^(q,j)  and  the  anharmonic  force 
constants.  The  errors  associated  with  the  u)^(q»j)  are  difficult  to  assess, 
but  it  can  be  stated  that  our  lattice-dynamical  calculations  are  in  close 
agreement  with  available  low-temperature  dispersion  curves,  reported  from 
inelastic  neutron  scattering  experiments,  with  the  most  extreme  differences 
being  less  than  107.  over  small  regions  of  the  Brillouin  zone,  and  that  an¬ 
harmonic  calculations  for  a  few  compounds  using  ua^(q,j)  values  differing 
uniformly  by  10%  throughout  the  Brillouin  zone  lead  to  only  107.  changes  in 
the  real  and  imaginary  anharmonic  self-energies .  The  anharmonic  force  con¬ 
stants  are  sensitive  to  the  input  values  of  tq  and  p,  and  in  both  cases 
values  of  these  at  low  temperatures  are  needed.  Values  of  0  are  typically 
determined  directly  at  4  K  with  an  error  of  about  ±  2%.  Such  an  error  in  B 
will  lead  to  errors  of  271  or  less  in  the  various  coefficients  associated 
with  the  anharmonic  force  constants.  Although  rQ  is  usually  determined 
quite  precisely  at  290  K  to  better  than  ±  0.17.  via  X-ray  measurements, 
values  of  rQ  at  lower  temperatures  are  often  not  directly  utaasured  but  are 

determined  from  the  known  290  K  value  of  r  and  the  known  thermal -expans ion 

o 

data.  This  leads  to  a  typical  error  of  about  1  0.27.  in  rQ  at  4  K.  such 
an  error  will  generate  errors  of  between  6-107.  for  the  various  coefficients 
associated  with  the  anharmonic  force  constants.  The  accumulated  errors 
in  these  coefficients  will,  of  course,  lead  to  errors  in  the  individual 

contributions  to  the  different  total  self-energy  quantities.  It  should  be 

A  (47  (87 

recalled  that  ApCOjt.fi)  is  obtained  as  a  sum  of  A  (0,t),  A  (0,t), 

A<6)(0,t,C7)  and  A^\o»t»Q),  and  although  the  errors  on  the  individual  con¬ 
tributions  to  the  self-energy  stemming  from  errors  in  rQ  and  S  are  no  larger 

A 

than  about  107.,  the  accumulated  error  on  A^(0,t,f7)  is  considerably  larger 


A 


because  of  the  opposite  signs  of  these  individual  contributions  and  because 

of  the  often  close  equivalence  of  the  two  dominant,  but  competing,  terms 

A^)(0,t)  and  A^)(0,t,fl).  However,  r(0,t,fi)  is  simply  the  sum  of  two 

components,  1*^(0, t,Q)  and  rW(0,t,fl),  both  of  which  are  accurate  to  about 

+62,  and  hence  the  accumulated  error  on  r(0,t,fi)  will  be  smaller  than  that 

for  AA(0,t,fl).  In  Tables  5-8  we  have  listed  our  estimates  of  the  possible 

errors  associated  with  the  different  calculated  self-energy  components  based 

on  the  known  errors  on  r  and  $  only. 

o 

E 

The  errors  on  the  experimental  determinations  of  A  (0,t), 

A 

AT<0,t,fi),  and  r(0,t,wt)  can  vary  significantly.  The  errors  on  the  measured 

w.j,  Q<0,t)  are  generally  quite  small  and  at  worst  no  larger  than  about  +1%, 

but  the  errors  on  w  (0,t)  will  be  somewhat  larger  than  this  due  to  the 
*  »  * 

error  in  calculating  P.  The  experimental  and  theoretical  results  strongly 
suggest  that  the  zero-temperature  anharmonic  self-energies  are  quite  small, 
so  that  the  use  of  the  approximation  w  oof0’c)  *  w^(0,t)  should  not  generate 
errors  of  more  than  a  few  percent  in  the  determinations.  Consequently  we 

believe  that  the  determinations  of  AE(C,t)  and  A' (T)  are  generally  accurate 

T  A 

to  better  than  +1.5  cm”  .  Of  course,  the  determinations  of  A^CO.t.w^)  are 

somewhat  worse  than  this  depending  on  the  accuracy  of  the  extrapolation 

procedure  used  to  determine  AA(0,t ,ui  ’).  The  errors  on  the  measured 

0  t 

r(0,t,w  )  are  determined  to  be  no  better  than  +15%  for  the  longei  phonon 
lifetimes  and  no  better  than  +25%  for  the  shorter  phonon  lifetimes. 

In  conclusion,  therefore,  we  have  found  a  good  agreement  between 
theory  and  experiment  concerning  the  sign,  magnitude  and  temperature  depend- 
ence  of  q  ■v  0  transverse  optic  phonon  self-energies  in  a  range  of  weakly 
anharmonic  systems.  Although  discrepancies  do  exist,  they  are  generally 
within  Che  errors  associated  with  either  the  calculations  and/or  the  exper- 


iments.  For  the  range  of  temperatures  and  frequencies  covered  by  our  work, 
therefore,  the  current  theories  of  weak  anharmonicity  appear  to  satisfac¬ 
torily  describe  the  anharmonic  interactions  in  materials  like  simple  cubic 


ionic  solids. 
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III.  STRONGLY  ANHARMONIC  SYSTEMS:  DISP1.ACIVE  FERROELECTRICS 

A.  Introduction 

Over  the  past  two  decades  the  existence  of  a  soft  mode  has  been  observed 
in  many  ferroelectric  materials.  Soft  mode  behavior,  that  is  the  rapid  move¬ 
ment  to  lower  frequencies  by  a  lattice  mode  as  the  temperature  is  decreased, 

(37)  (38) 

was  first  proposed  by  Cochran  and  Anderson  to  explain  the  anomalous  tempera- 

i 

ture  dependence  of  the  dielectric  response  of  ferroelectric  materials. 

As  discussed  in  the  earlier  sections,  the  temperature  dependence  of  any 
lattice  mode  frequency  arises  through  the  anharmonic  terms  in  the  lattice 
potential  energ  y  •  As  shown  in  the  previous  section,  the  hermitean  self-energy 
contributions  to  the  normal  mode  frequencies  in  most  systems  are  smaller  than 
the  harmonic  contributions.  In  displacive  ferroelectrics ,  however,  it  was 
theorized  that  the  anharmonic  contributions  to  the  soft  mode  were  larger  than 
the  harmonic  contributions  and  were  thereby  responsible  for  stabilizing  the 
paraelectric  states.  In  this  section  ve  describe  our  experimental  measure¬ 
ments  of  the  temperature  and  frequency  dependence  of  the  self-energies  asso¬ 
ciated  with  the  soft  mode  in  displacive  f erroelectrics  which  were  the  first  to 

directly  confirm  these  theoretical  predictions. 

A  TT 

B.  Experimental  Results  for  A  (o,t,0)  and  A"(o,t) 

( jc,) 

Cowley'  has  discussed  the  stability  of  an  anharmonic  crystal  using 
the  techniques  of  quantum  field  theory  and  has  analyzed  the  dielectric  response 
of  such  an  anharmonic  crystal.  The  lattice  contribution  to  the  complex  di- 


we  shall  assume  that  the  dielectric  response  is  determined  prime  ily  by  the 
summation  term  in  equation  (38)  associated  with  the  soft  mode  and  shall  ne¬ 
glect  contributions  from  the  other  q  c  o  transverse  optic  phonons  and  any 


temperature  dependence  of  the  f  .  The  self-energy  D(o,j,fl)  is  defined  by 

equations  5-18  of  the  previous  section. 

At  0  ■  o,  the  self-energy  terms  associated  with  the  soft-mode  can  be 

related  to  experimental  parameters  involving  the  lattice  contribution  to  the 

(AO) 

static  dielectric  response  as  follows:  ' 


2tth(o.t)AE(o,t) 


AeL(o)l 


a>b(o,t)2  +  2mhti,t)A^  (o,t,o)  AeL(o)l  +  (o) 

^  P  xi 


a)h(o,t)2  +  2(jh(o,t)A^(o,t,o)  tTi(°) 
u>h(o,t)2  +  2u)h(o,t)A^  fot.o)  (o)+Ae^(o)]v 


=  A(T) 


where  (o,t)  d -.notes  the  soft  mode,  Ae^(o)]^  is  the  change  in  eL(o)  in  raising 
the  temperature  from  to  TK  isobarically  at  1  bar,  AeL(o)]T  is  the  change  in 
eL(o)  in  raising  the  pressure  isotliermally  at  a  temperature  T  so  as  to  restore 
the  volume  of  the  crystal  to  that  which  it  has  at  T^K  and  1  bar  pressure,  and 
Ac^(o)]v  is  the  change  in  eL(o)  on  raising  the  temperature  to  If.  isobarically 
at  the  crystal  volume  at  T^.K  and  1  bar  pressure.  The  equations  show  that 
suitable  measurements  of  the  temperature  and  pressure  dependence  of  c*J(o) , 


therefore,  will  lead  to  assessments  o*  the  soft  mode  self-energy  components. 


Figures  15  and  16  show  the  results  of  such  measurements  for  the  tempera¬ 
ture  dependence  of  the  soft  mode  self-energy  components  for  KTaO^  and  SrTiO^ 
h  2 

An  estimate  of  ui  (o,t)  can  be  made  from  this  data.  At  sufficiently  high 
temperatures,  if  the  thermal  population  factors  contained  in  the  theoretical 
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f ^ 

expression  for  A  (o,t,o)  are  expanded  in  terms  of  tlm/kT  <  1,  it  is  found  that 

A 

the  leading  terms  in  both  the  cubic  and  quartic  contributions  to  A  (o,t,o)  are 
linear  in  T.  The  temperature  dependence  of  the  function  A(T)  is  shown  in 
figures  15  and  16  and  does  appear  to  be  reasonably  linear  with  tempera- 
ture  thereby  implying  a  linear  dependence  of  A  (o,t,o)  on  T.  An  extrapola¬ 
tion  of  this  linear  part  back  to  T  =  o  should  therefore  yield  an  interept  of 

approximately  w  (o,t)  /u  (o,t)  .  Using  measured  low-temperature  values  of 

*i 

u>_.  (o  *  t) 2  of  o'  (o ,  t) 2  =  0.35  X  1024  H„  for  KTaO,  and  a),  ..(o  ,t)2  •  1.79  X  1024H„ 
o  2  3  100  2 

for  SrTiO^  yields  values  of  to^(o,t)2  =  -0.004  X  1024  A^  for  KTaO^  and 
wh(o,t)2  =  -1.43  X  1024  for  SrTi03. 

Our  results  imply,  therefore,  that  the  soft  mode  harmonic  frequency  in 

SrTiO^,  and  probably  in  KTaO^,  is  imaginary  and,  furthermore,  is  stabilized 

h  A 

by  the  large  anharmonic  self-energy  components  2w  (o,t)A  (o,t,o)  which  domi¬ 
nate  all  the  self-energy  contributions  combined. 

A 

C.  Experimental  Results  for  A  (o,t,ft)  and  r(o.t.ft) 

Following  the  theoretical  and  experimental  analyses  presented  in  sec¬ 
tion  II,  we  have  also  determined  the  frequency  deperdence  of  the  self-energy 
terms  2w^(o,t)A^(o,t ,fi)  and  2w^(o,t)r(o,t,fi)  via  Dispersive  tourier  Transform 
Spectroscopy.  The  frequency  dependence  of  the  self-energies  determined  in 
this  way  are  shown  for  UTaO^  in  figures  17  and  18.  The  results  reveal 
again  that  these  soft  mode  self-energies  are  extremely  large  compared  to  those 
found  Jr.  weakly  anharmonic  systems.  In  addition,  these  self-energies  have  a 

strong  temperature  dependence  throughout  the  frequency  range  covered.  It  is 

h  A 

of  interest  also  to  note  that  the  determinations  of  2u>  (o,t)A  (o,t,ft)  at 
ft  «  o  determined  from  the  static  dielectric  constant  measurements  discussed 
in  part  B  of  this  section  are  closely  consistent  with  the  values  determined 


via  the  DFTS  measurements. 


IV.  HYDROGEN  BONDED  FERROELECTRICS 


A .  Intr oduct ion 


Recent  discussions  of  the  transition  in  hydrogen-bonded  ferro¬ 
electric  systems  have  examined  the  consequences  of  interactions  between  the 

ferroelectric  mode  and  fluctuations  in  the  phonon  density.  Coombs  and 
/*>\ 

Cowley have  used  a  weakly  anharmonic  theory,  within  the  framework  of  the 
Landau  theory  of  phase  transitions,  to  consider  the  coupling  of  the  ferro- 
electrJc  mode  to  fluctuations  in  the  acoustic  phonon  density,  and  they  find 
a  response  function  for  the  ferroelectric  mode  of 


Gi(uj) 


uujy^  -  u) 


[E  Li/(1 


-  iu)T>] 


Here  oj^  is  the  harmonic  ferroelectric  mode  wavenumber,  the  renormalized 
ferroelectric  mode  wavenumber  containing  conventional  anharmonic  self-energy 
contributions,  y^  the  mode  damping,  t  an  average  lifetime  for  the  acoustic 
phonon  density  and  is  given  by 


7Zu> 

Li  =  kFnv(0’1’~1)l  V"i  +  1) 


where  the  V  coefficients  are  second-order  cubic  anharmonic  coupling  coef- 

(4) 

ficients  and  the  are  thermal  population  factors.  Young  and  Elliott 
have  extended  the  pseudospin  model  of  Kobayashi^*^  to  include  coupling 
between  the  soft  pseudospin  mode  and  fluctuations  in  the  phonon  density  via 
two-phonon  terms  in  the  Hamiltonian,  and  they  find  a  response  function  for 
the  ferroelectric  mode  of 


G^tu) 


_ _ _ i _ 

h^(u)^  -  ituy^  “  is^)  -  [4r(aX)L/(l  -  iurr)] 
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Here  T  is  the  energy  of  the  pseudospin  mode,  (o  )  the  thermal  average  of 
the  spin  operator  and  L  is  given  by 


1  c’  2  2lMn,-  +  1> 

l’h  l  -ns —  • 


1,2 


(44) 


The  responses  from  these  two  different  theoretical  approaches  are  structurally 

identical,  but  differ  from  the  characteristic  response  of  normal  phonons  by 

the  addition  of  a  relaxing  self-energy  contribution.  This  has  two  important 

consequences:  it  predicts  the  existence  of  a  quasi-elastic  central  component 

and  the  divergence  of  the  ferroelectric  mode  susceptibility  at  a  temperature, 

T^,  below  the  characteristic  clamped  Curie-Weiss  temperature. 

Attempts  to  determine  the  magnitude  of  the  relaxing  self-energy 

have  been  made  for  a  number  of  hydrogen-bonded  ferroelectrics  via  experiments 

to  measure  T.  and  via  searches  for  the  central  component.  The  former  suffer 
A 

from  the  fact  that  the  ferroelectric  mode  wavenumber  is  determined  from 

extrapolations  of  data  which  are  deduced  from  fits  to  a  simplified,  and 

possibly  ambiguous,  model  involving  the  coupling  of  the  ferroelectric  mode 
(42-44) 

to  other  phonons.  Such  experiments  generally  have  led  to  the  con¬ 

clusion  that,  if  present,  the  relaxing  self-energy  is  quite  small,  but  the 
inaccuracy  of  the  method  often  precludes  a  definitive  statement  as  to  whether 

(44) 

the  relaxing  self-energy  is  finite  or  not.  Although  Lagakos  and  Cummins 

have  reported  on  a  preliminary  observation  of  a  central  component  in 

(KDP)  which  also  suggests  that  the  relaxing  self-energy  is  quite  small, 
(45) 

the  same  authors  have  found  no  evidence  for  a  central  component  in 
CsHjAsO^  (CsDA)  . 


« 
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The  possibility  of  a  separate  determination  of  the  magnitude  of 

(4) 

the  relaxing  self-energy  has  been  pointed  out  by  Young  and  Elliott.  They 
have  commented  that  the  inclusion  of  spin-two-phonon  coupling  in  their 
Hamiltonian  also  causes  a  splitting  of  the  E-modes  in  the  ferroelectric  phase 
where  the  corresponding  singly  degenerate  modes  have  symmetries  of  B^> 

This  splitting  will  be  in  addition  to  any  reststrahlen  splitting,  and  may 
be  distinguished  from  it  because  their  combined  effect  is  dependent  on  the 
propagation  direction  of  the  phonon.  They  predict  the  splitting  to  be  given 
by 

id  2  -  <r  2  =  S2  for  k  along  (0,0,1)  (45) 

A  D 

o)A2  -  iDg2  =  S2  +  R2  for  k  along  (1,1,0)  (46) 

-  uig2  =  (S^  +  R^)^  for  k  along  (1,0,0)  (47) 

2 

where  S  is  the  relaxing  self-energy  splitting, 

,  4m  <aZ>V(0,l,2) 

S2=-^— ^ -  ,  (48) 

2 

and  R  is  the  reststrahlen  splitting, 


2  2  2 
R  =  "  “r 


(49) 


where  and  ajj.  are  the  zone-centre  longitudinal  and  transverse  optic  wave¬ 
numbers,  respectively. 

Suitable  experimental  measurements  of  any  E-mode  splitting, 
therefore,  can  lead  to  separate  estimates  of  the  magnitude  of  the  relaxing 
self-energy.  Such  measurements  can  also  lead  to  important  information  con¬ 
cerning  the  central  component  because  the  anharmonic  coupling  coefficients 
that  determine  the  magnitude  of  the  E-mode  splitting  are  also  related  to  the 


behavior  of  the  amplitude  and  wavenumber  of  the  central  component. 
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The  following  sections  describe  laser-Raman  investigations  by  this 
group^^  of  such  E-mode  splitting  in  hydrogen-bonded  ferroelectrics ,  the 
results  of  which  suggest  that  the  relaxational  processes  responsible  for  the 
central  components  have  sufficiently  long  lifetimes  as  to  make  their  obser¬ 
vation  feasible  through  electronic,  rather  than  optical,  techniques.  The 
final  sections  describe  such  a  search  via  radio-frequency  dielectric  constant 
measurements . 

B.  E-Mode  Splitting  Investigations 

Figure  19  shows  schematically  the  laser-Raman  scattering  geumetries 

used  to  study  the  E-mode  splitting. 

Figure  20  shows  in  detail  the  low  low-frequency  B^,  B2  spectra 

* 

recorded  in  the  ferroelectric  phase  of  KDP,  KDA  and  KD  A  for  the  three 

scattering  geometries  consistent  with  equations  (45)-(47),  and  figure  21 

* 

shows  the  full  B^,  E2  spectra  for  RbDA  and  RbD  A  for  tvo  of  the  scattering 

geometries.  As  these  figures  show,  many,  but  not  all,  of  the  E-modes  do  show 

definitive  splittings  arising  from  non-reststrahlen  effects.  For  reasons  we 

have  described  elsewhere , these  splittings  cannot  arise  from  eithei 

strain  or  birefrigence  artifacts,  and  hence  we  believe  the  splittings  to 

originate  from  the  relaxing  self -energy. 

The  response  function  given  in  equation  (43)  leads  to  the  conclusion 

that  the  ferroelectric  mode  will  condense  out  at  a  temperature  T,  such  that 

A 

X  X 

T^  <  Tc  ,  where  T£  is  the  clamped  Curie  temperature.  A  measure  of  the 
relative  strength  of  the  relaxing  self-energy  to  the  mode  total  self-energy 
is  given  by  a ,  where 
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a 


(50) 


This  response  function  also  leads  to  a  Debye  relaxation  form  for  the  central 
component  of 


Im  G(iu) 
G(0) 


4r(oX)L  tut* 

.2-2  ,  .  2  ' 

tl  0)  1  +  0)  T 


(51) 


where 


* 

T 


*2-2- 

■  0)  T 


1? w"  -  4F(oX)L 


_1_ 

U) 

c 


and  <uc  is  the  central  component  wavenumber.  In  the  high-temperature  limit, 
equation  (51)  leads  to  the  conclusion  that  u)c  will  have  a  temperature 
dependence  given  by 


(1  -  a) (T  -  T  X)  1 

u> - ~ - £ - 

C  [T  -  TcX(l  -  c)]  ^ 


(52) 


whilst  the  amplitude  of  the  central  component  will  be  given  by 


where  A  *  l/kTc  and  €  is  given  by 


153) 


(54) 


Analogous  results  to  those  given  in  equations  (52)  and  (53)  can  also  be 

(3) 

obtained  from  the  response  function  given  by  Coombs  and  Cowley  for  the 

ferroelectric  mode.  However,  it  is  to  be  noted  that,  although  the  two 

x 

theories  both  require  the  ferroelectric  mode  to  condense  out  below  Tc  ,  the 


temperature  at  which  this  occurs  is  not  necessarily  the  same  because  of  the 
different  temperature  dependence  predicted  for  the  ferroelectric  mode  by  the 
two  theories. 


The  lack  of  information  on  t  prevents  any  utilization  of  equation 

(52)  to  assess  the  specific  temperature  dependence  of  u^.  However,  estimates 

of  the  temperature  dependence  of  the  central  component  amplitude  can  be  made 

x 

via  equation  (53)  providing  measurements  of  T£  and  o  are  available. 

T  X  is  readily  determined  with  good  precision  via  suitable  dielec- 

(4) 

trie  or  elastic  constant  measurements .  or  is  given  by  aquation  (50)  and  as 


a 


ym  o,i, -m2 
L  V  hu>.  y 


(55a) 


a 


(55b) 


where  6u)^  is  a  mode  splitting  from  the  spin-two -phonon  interactions  for  k 
along  (0,G,1)  as  given  in  equation  (45).  Although  a  can  be  determined  via 
equation  (50),  as  mentioned  earlier,  such  calculations  are  not  reliable 
because  of  the  inaccuracies  in  determining  T^.  The  fundamental  difficulty 
with  equations (55a)  and  (55b)  is  that  they  cannot  be  applied  to  those  modes 
which  are  not  doubly  degenerate  in  the  paraelectric  phase  and  for  which 
splittings  may  thus  not  be  observed.  This  requires,  ultimately,  that  a  be 
estimated  using  equations  (45)-(47)  and  (55)  and  a  means  of  attributing 
values  to  V(0,l,-1)  for  those  modes  where  measurement  is  not  possible. 

Table  12  lists  values  of  o’  calculated  in  a  number  of  different  ways,  is 
determined  using  only  the  measured  E-raode  splittings  in  equation  (55b).  a^. 

is  determined  via  equations  (48)  and  (55a),  by  using  the  most  consistent 
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2 

values  of  S  determined  from  equations  (46)-(47).  It  is  to  be  noted  that 

there  is  little  difference  between  and 

fits  of  the  experimental  data  to  the  form  of  equations  (45)-(47).  (a)  and 

(cy)^  are  calculated  on  the  assumption  that  all  48  phonon  branches  contribute 

equally  to  a;  (a)  uses  the  observed  srlittings  averaged  over  all  13  allowable 

E-modes,  whilst  (a)  uses  the  splittings  averaged  over  only  the  observed 
m 

split  modes. 

Table  12  also  contains  estimates  of  a  determined  via  equation  (50) 

using  data  determined  from  studies  of  the  temperature  dependence  of  the 

x  (^7) 

ferroelectric  mode.  Using  measured  values  of  Tc  (from  Brody  and  Cummins 

for  KDP  and  from  Spillman  and  Lowndes for  the  arsenate  family),  or  and 

& 

a  were  then  calculated  using  appropriate  values  of  T  determined  by  fitting 
P  A 

the  measured  temperature  dependence  of  the  ferroelectric  mode  of  these  mate- 

(43) 

rials  (Lowndes  et  al.,  '  J.  F.  Ryan  19/2  (unpublished)  quoted  by  Lagakos 

and  Cummins  ^T^)  to  the  temperature  dependence  predicted  by  the  theories  of 

Coombs  and  Cowley^^  and  Young  and  Elliott  respectively. 

Inspection  of  Table  12  reveals  a  number  of  differences  between  the 

two  basic  approaches  to  determining  'ctinates  of  a.  Firstly,  the  a  values 

determined  from  the  E-mode  splittings  aie  consistently  smaller  than  the  values 

determined  from  studies  of  the  temperature  dependence  of  the  ferroelectric 

mode.  Secondly,  for  a  given  anion,  the  a  values  determined  from  the  E-mode 

splittings  increase  with  increasing  mass  of  the  cation,  whereas  the  reverse 

is  true  for  a  and  or  .  Thirdly,  the  a  values  determined  from  the  E-mode 
a  p 

splittings  are  always  positive,  whereas  a  and  a  are  negative  for  the  deu- 

a  p 

terated  arsenates.  As  mentioned  earlier,  T^  is  difficult  to  determine 


°E 


and  this  reflects  the  close 


precisely,  and  this  is  especially  true  for  the  deuterated  arsenates  where  two 
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(43) 

or  more  low-frequency  optical  phonons  are  coupled  to  the  ferroelectric  mode, 
thus  casting  doubt  on  the  adequacy  of  the  two  coupled  oscillator  fit  used  to 
determine  the  ferroelectric  mode  spectra  characteristics  from  the  Raman 
spectra;  we  believe  this  is  the  reason  for  T.  >  T  and  hence  for  the  negative 

A  C 

values  determined  for  <y  and  of  for  the  deuterated  arsenates. 

a  p 

Because  of  these  and  other  difficulties  in  precisely  determining 


the  temperature  T  at  which  the  ferroelectric  mode  condenses  out,  and  because 

A 

T.  and  T  are  apparently  very  close  for  the  materials  considered  here,  we 
Ac  ' 

believe  that  the  a  and  a  values  provide  less  reliable  estimates  of  or  than 
a  p 

those  determined  from  the  E-mode  splittings.  Although  these  latter  do  not 


lead  directly  to  an  estimate  of  the  total  contributions  to  a,  they  do  provide 

a  basis  for  establishing  a  range  of  possible  values  for  a.  For  instance,  a 

b 

provides  a  reasonably  precise  lower  bound  on  a  since  actual  observed  contri¬ 
butions  ace  included  in  its  calculation,  .’bilst  (a)  provides  a  reasonable, 

m 

but  less  distinctive,  upper  limit  on  a. 


Figure  22  shows  the  results  of  calculating  the  temperature  depen¬ 
dence  of  the  central  component  amplitude  for  the  hydrogenated  materials  using 

<y_ ,  ( a )  and  (a)  for  a  in  equation  (53).  Figure  23  shows  the  results  of 
Em 

similar  calculations  for  the  deuterated  arsenates.  Although  the  quantitative 


details  vary  somewhat  depending  on  the  particular  value  of  a  used,  the  quali¬ 
tative  trends  are  very  similar  irrespective  of  the  or  value  used.  The  central 
feature  of  importance  to  emerge  from  these  figures  is  that,  although  the 

X 

central  component  amplitude  grows  strongly  as  r  is  approached,  in  all  cases 
this  growth  is  curtailed  by  the  onset  of  the  ferroelectric  transition.  The 
data  indicate  that,  at  the  characteristic  transition  temperature,  the  largest 
growth  in  the  central  component  amplitude  occurs  for  KDP  and  the  smallest  for 
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CsDA.  This  may  be  an  important  contributing  factor  in  explaining  why 
(44  45) 

Lagakos  and  Cummins  *  were  able  to  find  evidence  for  a  central  component 

in  KDP  but  not  in  CsDA.  In  commencing  on  their  failure  to  find  a  central 

(45) 

component  in  CsDA,  Lagakos  and  Cummins  have  derived  two  upper  limits  for 

a:  firstly,  a  value  of  10  for  an  unresolvably  narrow  central  component, 

-2 

and,  secondly,  a  value  of  10  for  a  central  component  linewidth  exceeding 
the  free  spectral  range  of  their  Fabry-Perot  interferometer.  The  values  of 
a  determined  here  from  the  E-mode  splittings  are  nicely  bracketed  by  these 
limits,  but  tend  to  support  the  larger  value. 

C .  Radio  frequency  Dielectric  Constant  Investigations 
a)  Experimental  Studies 

The  light  scattering  measurements  of  the  E-mode  splittings  in 
hydrogen-bonded  ferrcc lectrics  described  in  the  preceding  section  suggest 
that  a  central  component  may  exist  in  these  materials  but.  at  frequencies 
which  may  be  too  small  to  be  resolved  by  optical  measurements.  Recent 
acoustic  alce.ruation,^^  electron  spin  resonance and  electron  para¬ 
magnetic  resonance^^  work  also  suggest  that  such  re laxat ional  processes 
do  occur  but  with  lifetimes  sufficiently  long  as  to  make  their  observation 
feasible  through  electronic,  rather  than  optica],  techniques.  We  have 
therefore  searched  for  the  existence  uf  central  components  via  radio¬ 
frequency  dielectric  constant  measurements. 

In  the  paraelectric  phase  of  KDP,  the  dielectric  response  measured 
in  the  5-80  MHz  range  exhibits  no  frequency  dependence  which  would  be  charac¬ 
teristic  of  a  relaxation  process  and  generally  displays  a  temperature 
dependence  commensurate  with  measurements  at  lower  frequencies.  Within  the 
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sensltlvity  of  our  experiments,  therefore,  no  significant  evidence  was  found 
for  any  contribution  to  the  dielectric  response  in  this  frequency  range 

which  arises  from  the  existence  of  a  central  component. 

(52) 

However,  as  we  have  reported,  the  response  for  CsDA  in  the 
paraelectric  phase  does  exhibit  a  marked  frequency-dependent  dielectric  loss, 
especially  at  the  higher  frequencies,  for  temperatures  just  above  the  tran¬ 
sition.  Unfortunately,  the  characteristic  frequency  of  the  mechanism 
associated  with  the  dielectric  loss  appears  to  be  above  the  80  MHz  upper 
limit  or  our  applied  field  frequency  using  this  measuring  technique.  Although 
e7  shows  a  continuous  decrease  with  increasing  temperature  from  the  transi¬ 
tion  at  149  X  5K,  c" ,  especially  at  the  higher  frequencies,  first  grows 
rapidly  as  the  temperature  is  raised  from  the  transition  and  then  decreases 
with  a  further  increase  in  temperature  until  it  is  barely  measurable  at 
about  170K. 

It  is  difficult  to  comment  on  the  origins  of  the  mechanism  causing 
this  dielectric  loss  in  the  absence  of  a  full  scan  of  its  frequency  dependence. 
The  results  are  unlikely  to  stem  from  impurities  or  defects  in  the  sample 
which  typically  lead  to  relaxation  tines  several  orders  of  magnitude  longer 
than  those  that  could  be  involved  here.  Rather,  the  results  are  more  likely 
to  be  intrinsically  related  to  the  onset  of  the  transition,  in  view  of  the 
growth  of  this  dielectric  loss  as  the  transition  temperature  is  approached. 
Although  the  results  could  be  interpreted  as  being  consistent  with  the 
findings  of  Lagakos  and  Cummins  (in  the  sense  that  if  the  central  peak  for 
KDP  lies  at  a  higher  frequency  than  that  for  CsDA  then  the  optical  measure¬ 
ments  would  be  more  favorable  for  studying  the  central  peak  in  KDP  and 
the  dielectric  measurements  for  studying  it  in  CsDA),  the  present  data  for 
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CsDA  would  not  seem  to  be  described  by  a  simple  Debye  relaxation  form  with  an 
amplitude  growing  continuously  as  the  clamped  Curie  temperature  at  126  K  is 
approached,  as  is  suggested  by  the  current  theories. 

The  results  for  the  ferroelectric  phase  of  both  materials  also  reveal 
some  significant  dielectric  behaviour  which  has  not  previously  been  reported. 
For  both  materials,  the  dielectric  response  exhibits  a  pronounced  dependence 
on  frequency  which  is  well  described  by  a  Debye  form  of  relaxation: 

e( to)  *-  e'(o>)  +  (»)  =  e(“)  +  r— —r (56) 

i  1Cl)T 

where  e(w)  is  the  complex  dielectric  response  at  a  frequency  w,  e(°0  is  the 
complex  dielectric  response  at  u  ^100  MHz,  A  is  the  amplitude  and  t  the 
characteristic  relaxation  time.  Data  at  representative  temperatures,  to¬ 
gether  with  their  associated  fits  to  the  form  of  equation  (56),  are  shown  in 
figures  24  and  25.  Table  13  summarizes  values  of  the  parameters  characterizing 
the  Debye  relaxation  process  for  both  materials.  Although  the  dielectric 
relaxation  occurs  at  somewhat  lower  frequencies  for  K.rP  than  for  CsDA,  the 
relaxation  data  for  the  two  materials  are  qualitatively  similar  with  the  t 
decreasing  dramatically  with  increasing  temperature  in  the  range 

T  *  0.8  T  to  T  :  T  where  T  is  the  transition  temperature.  The  amcii- 
tr  tr  tr 

tude  also  increases  substantially  over  much  of  this  same  temperature 
range  but  appears  to  level  off  close  to  the  transition  temperature  for  KDP 
and  to  actually  go  through  a  maximum  about  5  K  below  the  transition  tempera¬ 
ture  for  CsDA.  Two  further  aspects  of  these  parameters  should  be  noted. 
Firstly,  e(»)  has  a  real  part  whose  temperature  depencence  is  weakly  sug¬ 
gestive  of  divergent  behaviour  but  which  is  not  described  by  a  Curie-Weiss 
behaviour,  and  an  imaginary  part  which  suggests  a  source  of  high-frequency 
conductivity  which  is  independent  of  the  relaxation  process.  Secondly, 
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values  of  t  are  well  fitted  by  the  thermal  activation  fora 

t  *=  t  exp(E/kT) 
o 

with  t  -  4.114  X  l<f10s  and  E  *=  0.0466  eV  for  KDP,  and  t  =  1.586  X  l(T10s 
o  o 

and  E  =  0.0603  eV  for  CsDA. 
b)  Theoretical  Discussion 

Almost  simultaneously  with  our  observation  of  the  central  component  in  the 

power  spectrum  of  homogeneous  polarization  fluctuations  in  the  ferroelectric  phase  of 

(53) 

KDP  and  CsDA  discussed  above ,  Memelstein  and  Cummins  reported  on  a  central  component 
observed  in  the  power  spectrum  of  polarisability  fluctuations  in  the  ferro¬ 
electric  phase  of  KDP. 

Our  subsequent  analysis  of  the  results  for  KDP  (and  by  implication  for 
CsDA,  etc.)  suggests  that  both  effects  may  be  attributed  to  a  homogeneous  expo¬ 
nential  decay  of  the  polarization.  In  our  dielectric  experiments,  the  response 
of  the  polarization  due  to  an  applied  electric  field, 

X(«>  -  <«7«£)u  +  iQ+  ,  (57) 

was  directly  measured  in  the  ordered  phase  of  KH^PO^.  In  equation  (44), 

6P  and  6E  are  in  the  direction  of  the  spontaneous  polarization.  As  reported, 
the  complex  dielectric  constant  was  well  fitted  by  the  simple  form 


X<»>  -  (£  -  ’  (T  <  V  (58a) 

in  the  (oj)  ~  10^  Hz  range.  The  fluctuation-dissipation  theorem  for  polariza¬ 
tion  fluctuations 


<i‘V 


v  yW'*1 

TT  U) 


58  b) 


may  therefore  be  rewritten  as  a  Lorentzian  central  component, 

lim  < 1 AP  | 2  -  (kRTxT/")[T/(l  +  (wt)2)],  (T  <  T  ). 
Q+o  Qw  —  B  T  c 


(59) 
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The  light-scatterin.  rate  in  the  low-frequency  regime,  u>  «  (k^T/fi),  is  well 
known  to  obey  „  3 

'  -  <^>  <l<J  2>  • 

where  w  =  c  |k^  |,  =  c|k^|,  and  e^,  represent  the  frequencies  and  polari¬ 

zations,  respectively,  of  the  initial  (i)  and  final  (f)  states. 


k.  -  k, ,  uj 

V  X  I 


the  polarization  matrix  elements  of  the  fluctuating  pol arisability  are 


**£  -5fx* 


X  e,. 


•^1 


(61) 


(62) 


and  c  is  the  speed  of  light  in  the  ferroelectric  medium.  Since  we  wish  to 
establish  the  extent  to  which  the  experimental  result  in  equation  09)  will  be 
observed  by  measuring  the  effect  given  by  equation (60),  we  need  to  consider  the 
total  light  scattering  into  a  given  solid  angle 

m .  4 


--  (?)  <l*x“  2>  • 


(63) 


With  F  as  the  free  energy  per  unit  volume  we  have 

dF  =  -SdT  +  E  ■  dP  +  ]*X  dx.,  (64) 

-  3 

where  Xj  and  are  the  stress  and  strain  components.  A  complete  assess¬ 
ment  of  the  conributior.s  to  light  scattering  arising  from  thermal  fluctuations 
will  require  the  following  thermal  response  functions  to  be  considered: 

(6  5a) 


V*  *  I<3S/3I)P,x  • 


“j  =  <8P/3Vt,E,X  -  <ay3E)T,E,X* 


and 


V  -  OyiVr.E.!  ■  <3i'k/SXj)I,E,X- 
The  following  are  valid  as  Q  -*■  0: 


(65  b) 
(65  c) 

(65d) 


■50 


■  'A,,1' 

(66a) 

A/  A/ 

<ATq aSq>  -  °» 

(66b) 

v 

<AT*Ax^>  =  0, 

<44  v  -  v*t  • 

(66  c) 

(66d) 

MS 

<AVxjQ>  "  kBT2j  » 

(66e) 

<AxjQAxkQ>  ■  kBTXjk  ■ 

(66f) 

Equations  (66)  represent  a  complete  list  of  thermal  fluctuations  which  can  con¬ 
tribute  to  light  scattering. 

Expanding  x*1  in  the  local  thermal  parameters  T,P  and  strains  {x.}, 

~  J 

we  get 

dxfl  -  r^dT  +  rj1  •  dP  +  [r^dx.  (67) 


which  defines  the  V  coefficients.  The  full  static  form  factor  from  equations 
(66)  and  (c.7)  is  therefore 


lim 
Q  0 


’(kBTJ/CP.x)  +  (l“>* 


’  ^Xp  ^ 


+ 


.rfi„fi.«. 

(Xp  ‘j  5  1 


k,Ta. 

h 


(68) 


which  may  be  substituted  in  equation  (63)  to  give  the  total  light-scattering 
cross-section.  To  generalize  the  result  to  frequency  power  spectra  via  the 
transformation 


I  fi|2  I  fi|2 

"IXq  I  >  -  "IXqJ  >  * 

A/  A/ 

one  must  replace  the  fluctuation  matrix  (equations  66) 


“■  “VV  '  kBTRIk’ 

Q  -*■  0  Z 


(69) 


(70) 


by  dynamic  fluctuations 


<AY  n  AY.  „ 
jQu  kQto 

A/  V 


[Rjk(Q,a>)  -  ^(Q.w)] 


leaving  the  T  coefficients  in  equation  (68)  unchanged.  The  dynamical  equiva¬ 
lent  of  equations  (65)  must  then  be  obtained: 

«£•“>  -  +  l0+  •  <723> 

*<!>•“>  •  (5I/sPq,„  +  io+  •  (72b> 

a  (Q»u>)  «  («x,/SE)ft  .  .  ,  (72c) 

~  j  Q,w  +  io' 

VS-“>  ’  “V'Vq,.  +  io+  ■  <72d» 

Since  the  derivation  of  the  ten  linearised  equations  required  to  evaluate 
equations  (72) is  not  necessary  here,  we  will  limit  our  evaluation  of  the  con¬ 
tributions  to  equation  (68)  to  an  estimation  employing  only  thermodynamic 
considerations. 

(53) 

Mermelstein  and  Cummins  attribute  the  central  component  observed 

in  their  measurements  to  a  thermal  diffusion  mode.  If  one  takes  only  the 

first  term  in  cur  equation  (68),  their  equation  (11)  is  the  result.  Hence 

they  effectively  write 

2.  2 
(~  )  (centra*-  component)  -  j  | 

•  (kBT2/™)IM(l/C(Q,u3)),  (73) 

with  a  dynamic  heat  capacity  of  the  thermal  diffusive  form 

[c  /C(Q,0J)  ]  -  [iD  q2/(«  +  iD_Q2)].  (74) 

rX  T~  T  ** 

Here,  ■  K/Cpx  is  the  thermal  diffusivity,  and  K.  is  the  heat  conductivity. 
This  yields  a  Lorenztian  central  component  of  width 

1/tq  *  dtq2  • 


(75) 


.  •  ■  mfr 
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I£  the  central  component  stems  from  the  second  term  of  equation  (68) 
we  may  write 

3 

k  T  r ,  ^  £. 

(central  component)  -  ( — r^)  (-7-)  ImKr*  )  •  x(Q»“>)  •  (Ik,1)],  (76) 

**  7TW  vr  —  ~  r 

C 


2  <-nrMrfl 


\V 

A/ 


(77) 


From  equations  (59)  and  0^7)  one  would  then  expect  a  Lorentzian  central  component 
with  a  Q  independent  width  of  1/t.  This  would  allow  the  observations  of  both 
experiments  to  be  described  by  the  same  microscopic  physical  process  which, 
however,  cannot  be  unambiguously  identified  at  this  time. 

First  l«t  us  consider  the  width  of  the  central  component.  The  light- 

2 

scattering  width  agrees  with  both  D^Q  and  the  dielectrically-measured  width 
(which  is  obviously  a  Q  -  0  homogeneous  mode).  One  has  two  possibilities: 

(i)  The  light-scattering  central  component  arises  from  a  thermal  diffusion 
mode  with  a  width  given  by  equation  (75),  while  the  polarisation  central 
component  arises  from  some  other  mechanism.  The  similarity  of  the  two  line- 
widths  would  then  be  pui  :ly  coincidental.  (ii)  Both  central  components  have 
the  same  physical  source.  Then  the  agreement  with  equation  (75)  is  purely 
coincidental.  This  dichotomy  could  be  resolved  experimentally,  of  course,  by 
changing  Q  in  the  light-scattering  work  and  observing  whether  or  not  the 
width  appreciably  changes. 

Secondly,  let  us  consider  the  relative  intensities  of  the  scattered 
light  given  by  the  first  and  second  terms  in  equations  (77).  The.-,*  arc 


h(thermal  diffusion)  «  (7— )|t^  •  -jr 

~Px  * 


w 


(78) 


2  it  It  £  It 

h(polarisation)  «  Jp(cf  *  XT  '  Cj)  •  ^p(ef  •  Xy  '  s^).  (79) 
An  analysis  of  equation  (78)  has  been  carried  out  by  Mermelstein  and 
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(53) 

Cummins  .  An  equivalent  analysis  of  equation  (79!  is  much  more  diffi¬ 
cult  in  that  (3/3P)x^,  is  a  third-order  tensor  which  has  not  (to  our  knowledge) 
been  directly  measured.  The  depolarization  effects  will  be  substantial  and 
not  qualitatively  different  from  equation  (77).  For  an  estimate  of  absolute 
size  we  take  the  Landau  approach  and  expand  the  free  energy  in  powers  of  P, 

(80) 


F(P,T)  =  l  a  (T)P 
n=l 


2n 


stopping  at  n  =  2  for  second-order  and  n  =  3  for  first-order  phase  transf. 

tions.  The  a  (T)  which  change  sign  at  T  *»  T  are  well  known^"^.  The  co- 
n  c 

efficients  are  fit  to  E  =  0  and  (3x^,/3P)  is  estimated  via  equations  (64)  and 
(80).  Although  this  procedure  is  only  a  tough  approximation,  it  allows  an 
order-of-magnitude  comparison  between  the  light-scattering  intensities 
given  by  terms  1  and  2  of  equation  (68) 


h(thermal  diffusion)  XT  - 
H(polarization)  3 


-2  -1 
10  to  10 


(81) 


This  suggests  that  the  homogeneous  polarization  effects  would  dominate  those 
of  a  thermal  diffusion  mode  if  it  were  present,  and  hence  supports  our  view 
that  the  central  components  observed  by  the  two  techniques  have  the  same 


origin. 
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The  high-pressure  far  infrared  cell  for  transmission  spectroscopy- 
in  the  frequency  range  3-400  cm"1  and  with  hydrostatic  pressures 
up  to  S  kbar.  A,  retaining  closure;  B,  window  mount;  C,  window; 
D,  packing;  E,  window  retraining  ring;  F,  gas  input  line. 


Fig.  2.  The  experimental  layout  for  high-pressure  far  infrared  spectroscopy. 
A,  Michels on  interferometric  module;  B,  high-pressure  cell  module; 

C,  detector  module;  D,  source;  F,  beam  splitter;  G,  moving  mirror; 

H,  high-pressure  cell;  I,  module  vacuum  window;  J,  detector; 

0,  mirrors. 


AMPLITUDE-PHASE  FOURIER  SPECTROSCOPY 


Fig. 


Fig.  4 


Fig.  3 


Fig.  4 


Experimental  layout  of  the  dispersive  Fourier  transform 
spectrometer. 


Schematic  view  of  the  sample  and  screens  used  with  the  dispersive 
Fourier  transform  spectrometer. 
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The  calculated  frequency  dependence  of  (a)  A  (0,t,Q)  and 
(b)  r(0,t,n)  for  NaBr  at  5  K  (•  •  •),  200  K  (-  •  -) , 

400  K  (-  ••  -),  600  K  (-  •••-),  800  K  (-  -  -)  and  1000  K  ( 


NaC!  105K 


Fig.  8.  The  experimental  (•  •  •)  and  calculated  ( - )  frequency 

dependence  of  AA(0,t,D)  for  NaCZ  at  105  K.  Also  shown 
are  values  of  AA(0,t,fi)  determined  at  Q  =  0  (0)  from  low 
frequency  dielectric  constant  measurements,  and  at  Q  *  yj 
from  far  infrared  transmission  measurements. 
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Fig.  9.  The  experimental  (•  •  •)  and  calculated  ( _ )  frequency 

dependence  of  AA<0,t,n)  for  RbCZ  at  300  K.  Also  shown  are 
values  of  AA(0,t,0)  determined  at  Q  =  0  (0)  from  low- 
frequency  dielectric  constant  measurements,  and  at 
*  (Ut  (A)  from  far  infrared  transmission  measurements. 


Fig.  10,  The  experimental  (•  •  •)  and  calculated  ( - )  frequency 

dependence  of  AA(0,t,Cl)  for  RbCi  at  105  K.  Also  shown 
are  values  of  A^(0,t,fi)  determined  at  fi  a  0  (0)  from  low- 
frequency  dielectric  constant  measurements,  and  at 
0  »  u)t  (A)  from  far  infrared  transmission  measurements. 
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Fig.  11.  The  experimental  (•  •  •)  and  calculated  ( _ )  frequency 

dependence  of  r(0,t,n)  for  NaCi  at  300  K.  Also  shown  are 
values  of  r(0,t,Q)  determined  at  n  *  id  (A)  from  far 
infrared  transmission  experiments. 
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Fig.  12.  The  experimental  (•  •  •)  and  calculated  ( - )  frequency 

dependence  of  r(0,t,n)  for  NaCZ  at  105  K.  Also  shown  are 
values  of  r(0,t,n)  determined  at  Q  -  u-  (A)  from  far 
infrared  transmission  experiments. 


Fig.  13.  The  experimental  (•  •  «)  and  calculated  ( _ )  frequency 

dependence  of  r(0,t,fj)  for  RbCi  at  300  K.  Also  shown  are 
values  of  r(0,t,Q)  determined  at  fl  -  (A)  from  far 

infrared  transmission  experiirents. 
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Fig.  14.  The  experimental  (•  •  •)  and  calculated  ( _ )  frequency 

dependence  of  r(0,t,fi)  for  RbCi  at  105  K.  Also  shown  arc 
values  of  r(0,t,0)  determined  at  fl  ■  (A)  from  far 

infrared  transmission  experiments. 


Fig. 


h  E 

16.  The  temperature  dependence  of  2a'  (0,t)  a  (0,t) 
and  [o^(0,t)2  +  2u£(0,t)A£(0,t  ,0)]  expressed 
a*  a  fraction  of  the  quasi-normal  soft  mode 
frequency  (squared)  at  0  K  for  SrTiC>3.  The 
broken  line  illustrates  the  extrapolation 
used  to  determine  o;h(0,t)2. 
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Fig.  20.  The  low-frequency  B]_,  B2  spectra  recorded  at 
10  K  for  phonons  propagating  in  the  (0,0,1), 
(1,0,0)  and  (1,1,0)  directions  for  (a)  KDP; 

■jt 

(b)  KDA ;  and  (c)  KD  A. 
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Fig.  23(a) 


Fig.  23(b) 


Fig.  23(c) 


Fig.  23.  The  temperature  dependence  of  the  amplitude  of  the 

central  component  for  KD*A,  RbD*A  and  CsD*A  calculated 
from  equation  (53)  using  (a)  or  =  q£,  (b)  a  *  (a),  and 
(c)  a  =  (or)  .  The  broken  lines  indicate  the  charac¬ 
teristic  tA 
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Fig.  24(a)  Fig.  24(b) 


Fig.  24.  The  radio frequency  dependence  of  (a)  e'  and  (b)  t " 

for  KDP  at  different  temperatures  in  the  ferroelectric 
phase.  The  full  curves  represent  fits  of  the  experi¬ 
mental  data  to  the  Debye  form  of  relaxation  given  in 
equation  (43).  A,  91.12  K;  •,  90.47  K;  o,  103.36  K: 
x,  110.46  K;  k,  122.25  K. 
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Fig.  25(a) 


Fig.  25(b) 


Fig.  25.  The  radiofrequency  dependence  of  (a)  e*  and  (b)  e" 

for  CsDA  at  different  temperatures  in  the  ferroelectric 
phase.  The  full  curves  represent  fits  of  the  experi¬ 
mental  data  to  the  Debye  form  of  relaxation  given  in 
equation  (43).  A,  123.85  K;  •,  128.76  K;  O,  133.78  K; 
x,  138.75  K;  i,  144.42  K. 
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Table  9«  Experimental  values,  in  cm*4,  ofA*(0,/)  and  AjkO,  determined  at  O  =0 
and  Q  «  a>,,  for  the  cesium  halides. 


A'iO.t) 

* 

TOO 

n-o 

Q  BU| 

A^O.t.O) 

AftO./.u,) 

s 

0 

A 

-1.1 

-0.7 

CsCl 

100 

-1.6 

-2.0 

1.2 

-0.1 

200 

-2.9 

-5.8 

3.5 

1.5 

300 

-3.8±0.4 

—  9.8  ±  1.0 

5.1  ±2.4 

2.3±  2.0 

s 

0 

0 

0.4 

-1.2 

CsBr 

100 

-1.2 

-1.6 

0.8 

-0.2 

200 

-3.2 

-4.2 

1.7 

0.6 

300 

— 5.2±0.6 

— 6.8±  0.7 

2.6±  1.5 

0.9±  1.5 

s 

0 

0 

0.2 

0.5 

Car 

100 

-1.0 

-1.4 

0.5 

1.3 

200 

-2.3 

-3.4 

1.0 

1.8 

300 

-3.7  ±0.4 

— 5.5±  0.S 

1.4  ±1.3 

2.2±1.1 
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Table  10. 

Experiment?.!  values,  in  cm' 

*,  otA*(0, 

')  and  ijtO.t.u), 

determined  at  fl 

and  O  »U|,  for  the  silver  and  thallium  halides. 

A*<0,t) 

T(K) 

0  =  0 

A^O.t.O) 

Afr0,t,u;t) 

5 

0 

0 

3.5 

-  4.5 

.rtgCl 

100 

-  0.8 

-1.7 

Q.3 

-  5.3 

200 

-  3.7 

-5.2 

-2.1 

—  7.2 

300 

-  8.8=  0.8 

-9.5*  1.0 

-3.0*  2.0 

-10.5*5.0 

5 

0 

0 

3.0 

-  1.4 

AgBr 

100 

-  0.6 

-1.3 

-0.3 

-  2.0 

200 

-  5.2 

-4.6 

-0.7 

-  3.2 

300 

-11.2*0.8 

-8.6  *0.8 

-1.2*  2.0 

-  4.8*  2.3 

5 

0 

0 

1.0 

3.0 

na 

100 

-  0.6 

—0.9 

3.9 

4.5 

200 

-  3.4 

-4.6 

7.1 

8.4 

300 

-  4.6*  0.7 

-8.3*  1.0 

10.8*2.0 

12.6  *2.1 

5 

0 

0 

1.0 

3.0 

100 

-  0.6 

-0.6 

3.1 

4.4 

200 

-  2.0 

-3.5 

6.0 

8.6 

300 

-  4.3*  0.6 

-7.6  ±0.9 

8.S*  2.0 

12.8*1.8 

TlBr 


Table  11.  Measured  values  of  the  temperature  and  pressure 

dependence  of  2r(0,t,u)t)  for  the  alkali  and  heavy 
metal  halides. 


T  (K) 

2 

2r(0.1,<()t)  (cm-1) 

90  200 

290 

<f2r(o.i  <j.)  , 

- L  (cm  1  kbar  ') 

290 

L1F 

7.3  *3.2 

10.4*2.9 

13.9*2.3 

17.1*1.9 

.  L1C1 

14.4*1.4 

17.4*1.4 

23.1*1.4 

34.5*1.4 

LiBr 

15.3*3.8 

20.5*4.7 

26.8*4.6 

34.6*5.4 

Lit 

25.0*4.7 

27.9*4.7 

35.0*7.5 

39.7  *8.8 

KaF 

6.S*1.3 

10.4*2.1 

12.7*2.1 

16.8*2.0 

0.10*0.07 

NiCl 

2.5  *0.7 

3.~±0.7 

5.3  *0.5 

7.0*0. 5 

-0.50*0.13 

NaBr 

3.4*0  9 

4.3  *0.7 

7.0  ±0.7 

8.4  *0.7 

-0.S2±0J4 

NaJ 

4.3  ±1.* 

6.1  *1.9 

9. 5*3.0 

13.7*2.9 

KF 

11.9*2.1 

14.0*2.1 

16.6*2.1 

19.0*2.4 

KCl 

1.4  ±0.6 

1.9  *0.4 

3.4  *0.5 

4. 5*0.4 

-0.08*0.07 

KBr 

0.9  ±0.5 

1.6  ±0.4 

2.8  *0.2 

4.7  ±0.2 

-0.03*  0.03 

KI 

2.0  *0.6 

27*0.5 

5.0  *0.5 

6.2  *0.5 

0.03  *  0.02 

RbF 

6.4  *0.9 

7.4  ±1.2 

9. 6*1. 3 

10.9*1.3 

RbCl 

o 

•II 

in 

cvi 

3.1  *0.4 

4.0  ±0.4 

4. 8*0. 4 

-0.04*0.03 

BbBr 

1.3*0. 3 

1.7  ±0.3 

2.4  *0.3 

3.0  ±0.3 

0.14*0.11 

Rbl 

1.1  ±0.2 

1.4  ±0.2 

2.1  *0.2 

2.8  *0.2 

0.12*0.04 

CsF 

9.2  *1.4 

9.9  ±2.0 

16.2*2.0 

19.7*2.0 

CsCT 

3.4*0. 5 

3.7  ±0.5 

5.1  ±0.5 

6.1  ±0.5 

-0.02*0.02 

CsBr 

1.5  =  0.2 

2.3  ±0.2 

3.1  *0.2 

4. 0*0. 2 

-0.04  *  0.02 

Csl 

1.0  =0.2 

1.4*0. 3 

2.1  *0.3 

2.6  *0.3 

-0.04*0.02 

AgCl 

3.5  =  0. 5 

6.0*0. 5 

9.5*0. 5 

13.2*1.0 

0.02  *  0.01 

A?Br 

2.5  =  0. 5 

4.5  *0.5 

8. 5*0.5 

13.2*  1.0 

0.04  *  0.02 

TIC1 

4.0  =  0.5 

5.0*0. 5 

6.3  *  0.5 

7.7  *  0.5 

-0.05*0.04 

HBr 

2.7  =  0.3 

3.0  *0.3 

3.5  *  0.4 

4.0*0. 4 

-0.03*0.02 
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Table  12.  The  different  values  of  a  determined  as  defined  in 
the  text.  Note  that  the  errors  quoted  for  as  and 
Op  reflect  only  the  inaccuracies  involved  in  the 
linear  extrapolation  of  the  ferroelectric  mode 
spectral  characteristics  and  do  not  contain  any 
estimate  of  the  inaccuracy  of  the  fitting  procedure. 


a,  x  101 

ap  x  10* 

c^x  10* 

oi  *  101 

<5>  x  101  <»>„  x  10! 

KDP 

6-7 

6-7 

020  ±  C02 

025  ±  003 

051 

206 

K.DA 

54  ±04 

306  ±31 

010  ±  002 

O10  ±  002 

030 

061 

RbDA 

3-6  ±02 

11  1  ±09 

025  ±005 

026  ±  005 

084 

1-88 

CsDA 

024  ±  OOl 

M  ±  01 

047  ±003 

052  ±  003 

HI 

251 

KD*A 

-5-6  ±04 

-68  ±  03 

039  ±  003 

043  ±  002 

070 

1-95 

RbD'A 

— 7-0  ±02 

-50  ±02 

068  ±  002 

075  ±  003 

1-50 

4-18 

CsD'A 

-7-3  ±01 

—5-4  ±  Ot 

065  ±002 

0  73  ±  003 

075 

301 
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Table  13.  Values  of  the  parameters  characterizing 

the  Debye  relaxation  process  in  the  ferro¬ 
electric  phase  of  KDP  and  CsDA  determined 
from  fits  of  the  experimental  data  to  the 
form  of  equation  (43).  The  transition 
occurs  at  122.4  K  and  149.5  K  in  KDP  and 
CsDA,  respectively 


T 

<K) 

«'(«) 

€"(ac) 

A 

1/t 

(MHz) 

122  35 

446-3 

386 

9146  ' 

3078 

12053 

3399 

51-0 

9174 

2789 

US  89 

283-9 

596 

895  3 

26-71 

117-13 

255-6 

60-1 

8924 

23-23 

115  55 

2791 

101-5 

796  7 

21  31 

KDP 

11446 

251-2 

927 

80S  1 

2073 

11046 

2308 

101-2 

763  1 

1807 

10693 

197-2 

882 

751-7 

15-63 

10336 

172-8 

952 

6178 

1339 

9847 

117-7 

452 

3400 

19  13 

91  12 

46-1 

165 

489 

18  84 

14947 

898 

17 

910 

75-39 

146-73 

397 

29 

1120 

52-60 

144  42 

313 

1-2 

122-7 

4504 

141-74 

25  7 

17 

121  1 

4152 

CsDA 

13875 

23-3 

4-8 

1115 

38C4 

133-78 

21-9 

7-1 

985 

3192 

128-76 

26-2 

109 

560 

2621 

126-41 

24-3 

9-7 

378 

24-78 

123  85 

21-9 

82 

166 

25-93 
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